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BOUNDARY VALUE AND EXPANSION PROBLEMS: ALGEBRAIC 
BASIS OF THE THEORY.* 


By R. D. CARMICHAEL. 


1. Introduction.—Linear problems in various forms have been central 
in the development of extensive chapters of mathematical analysis. Many 
of the profound phenomena of nature are subject to laws whose expression 
in mathematical form gives rise to fundamental linear problems of several 
kinds. Logically the simplest and historically the first to be treated in 
detail of the linear problems of pure mathematics are those having to do 
with systems of linear algebraic equations. And these hold the place of 
greatest importance both on account of their simplicity and complete 
development and on account of their suggestiveness in other linear problems. 
In fact it is true that a wide and important range of transcendental linear 
problems emerge from a direct consideration of the natural limiting cases 
of algebraic systems under the guidance of current problems of transcen- 
dental analysis. 

The deep-lying connection between algebraic and transcendental prob- 
lems of linear character has often afforded a useful help in the investigation 
of the latter, particularly in the case of integral equations. The object of 
the present investigation (of which the first portion is here presented) is to 
treat a considerable range of transcendental boundary value and expansion 
problems in close connection with certain algebraic problems of which they 
are the limiting cases. 

It is necessary first of all to formulate and solve the algebraic problems 
with special reference to their use as a heuristic guide in discovering truths 
and their proofs for the limiting cases which are matters of prime interest. 
From the outset classic theorems for certain transcendental cases throw a 
light back upon the algebraic questions and are of great value in formulating 
and in solving the algebraic problems. Thus we have from the outset a 
valuable interaction between the algebraic problems on the one hand and 
the transcfidental problems on the other. The great intimacy of this 
interaction becomes increasingly apparent as the investigation proceeds. 

In his researches on integral equations, integro-differential equations, 
functions of lines, and permutable functions, Volterrat has made frequent 


* Presented to the American Mathematical Society, April, 1920. 

t See his ‘‘Le cours sur les équations intégrales et les équations intégro-différentielles,” 
1913, where an exposition of his work is given with references to his earlier memoirs. 
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and consistent use of the guide to the transcendental problems afforded by 
looking upon them as limiting cases of algebraic problems. This seems 
to be the first instance of the systematic and profound use of the principle 
of guidance which in a modified or extended form is made to lead the way 
in our study of boundary value and expansion problems. The same 
general principle served Fredholm* as a heuristic guide in his fundamental 
memoir on integral equations; and Hilbertf carried through the limiting 
processes which served Fredholm merely as a guide and deduced his results 
as limiting cases of algebraic propositions when the number of variables 
becomes infinite. 

Though the guide to transcendental problems afforded by the properties 
of algebraic equations has attracted the greatest attention in the field of 
integral equations and the closely related matters, it was not here that the 
method was first employed heuristically. It is indeed the unpublished 
method which Sturm employed in his fundamental researches on differential 
equations of the second order. The differential equation with boundary 
conditions may be looked upon as the limiting case of a corresponding 
difference equation with boundary conditions, the latter being but an 
abbreviated form of a certain restricted system of algebraic equations. 
The passage to the limit from the difference equation to the differential 
equation may be carried through rigorously in the case of several boundary 
value problems, as was shown by Porter in 1902 (more than two years before 
Hilbert took a similar step for integral equations). 

An examination of the previous study of differential equations fromthe 
point of view in consideration brings out the fact that a given differential 
equation with boundary conditions may be realized in an infinite number 
of ways as a limiting form of algebraic systems. The investigator must 
therefore exercise care in choosing that algebraic system of which a maximum 
number of properties persist if he is to obtain the greatest advantage 
from the interaction of the two problems. ‘The results in § 6 of the present 
paper were suggested by classic theorems for differential equations; success 
in arriving at them required first the choice of a particular form of algebraic 
system to yield the differential equation and the generalization of this form 
so as to maintain the requisite characteristics. The scaffolding by which 
the results were discovered has not been reproduced in the section. 

Similarly, one may look upon partial differential equations, integro- 
* Acta Mathematica, 27 (1903): 365. 

Gottingen Nachrichten, 1904, p. 49. 

t For references and a fuller discussion of matters referred to in this paragraph the 
reader may consult Bécher’s lecture before the Fifth International Congress of Mathe- 
maticians and his “Les méthodes de Sturm.” 
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differential equations, and various types of mixed equations, each as a 
limiting form of appropriate algebraic systems. Moreover, one may con- 
sider sets of algebraic systems and their limiting forms of such sort as to 
furnish guidance in the study of systems consisting of differential equations, 
difference equations, integral equations, and so forth, and any one of many 
combinations of these. The algebraic theory developed in the present 
paper is of such sort as to furnish guidance in the investigation of boundary 
value and expansion problems associated with each of these various types 
of equations and systems, as I shall make clear in later papers recording 
the results of the present general investigation. 

The matter developed in this paper may be outlined as follows: In § 2 
are recorded preliminary considerations relative to sets of adjoint algebraic 
systems involving a finite number of parameters subject to determination 
so as to render the systems consistent. The conjugate character of the 
solutions of these systems is treated in § 3 and properties are established 
similar to those of orthogonality and biorthogonality in the theory of func- 
tions. These results are applied in $4 in deriving the ‘expansions’ of 
sets of constants in terms of conjugate sets. A certain graphic representa- 
tion of solutions is defined and analyzed in $5 and is employed in § 6 in 
deriving and extending certain theorems analogous to the Sturmian zero 
separation theorems for linear differential equations of the second order. 
In $7 there is a discussion of the general character of a fundamental 
limiting operation (which is to be employed frequently in analyzing the 
transcendental problems which are the goal of the investigation); and in 
§ 8 a second type of limiting case is briefly treated. Variation of parameters 
and Green’s functions for algebraic systems are the topics discussed in § 9. 
Finally, in $10, a method is given for the condensation of the algebraic 
expansions of § 4 into contour integrals, the methods and the results being 
analogous to those which Birkhoff (see reference in § 10) has developed for 
a certain class of differential systems. 

2. Adjoint Algebraic Systems.—Among the natural ways of introducing 
the so-called adjoint algebraic system in connection with a given algebraic 
system is one which is intimately related to the usual method of introducing 
the adjoint differential expression in connection with a given differential 
expression. We recall that the differential expression M(v) which is adjoint 
to a given homogeneous linear differential expression L(w) is gotten by 
seeking a v such that an indefinite integral of v£(u) can be written as a 
homogeneous linear differential expression in u of order one lower than 
that of L(u). Following up the analogy and employing a sum as to 2 in 
place of the integral used in the case of differential expressions, we are led to 
seek the conditions on y1, yz, «++, Yn such that for given a;; the coefficients 
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6;; exist of such sort that we have 


n n n 
i=1 j=l j=l 


for each k of the set 1, 2, ---, n. We see that we must have 


= 0 


¢=1 
for every k. Thus with the first of the two sets of linear forms 


n n 
j=1 


we are led to associate the second. If we start with the latter it is clear 
that we shall be led in the same way to the former. The two systems of 
forms are said to be adjoint each to the other. 
As an analogue of the classic Lagrange identity [involving vL(w) 

— uM(v) | in the theory of adjoint differential expressions, we have the 
readily verified identity 

n n n 

i=l j=l j=1 
which will be found to serve us in a fundamental way as a guide to the 


method of approach to numerous problems. 
Let us now consider the r homogeneous linear algebraic systems in the 


unknown quantities 2, 
nh 

(1) + + + =O, T= 1,2, +++, 


a separate system being formed for each value h of the set 1, 2, ---, 7. 
Here Aj, Ao, -- +, A, are 7 parameters by the choice of which the simultaneous 
consistency of the r systems is to be secured. With these systems let us 
associate the adjoint systems 


np 


j=! 


A necessary and sufficient condition that each of the r systems in (1) 
shall have solutions not identically zero is that each of the equations 


(3) | + + = 0, h=1,2,---,7, 


shall be satisfied, where for a given value of hf the first number is a deter- 
minant of order n;, whose element in 7th row and jth column is that which 
is written out explicitly. Obviously we have the same necessary and 


| 
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sufficient condition that each of the r systems in (2) shall have solutions 
not identically zero. 

A set of (finite) values of \i, Xe, «++, Ay for which (1) [(2) ] has for each 
h a solution not identically zero will be called a set of characteristic values 
for (1) [(2)]. We have just seen that (1) and (2) have the same sets of 
characteristic values, since the sets of characteristic values are the (finite) 
solutions of system (3) considered as a system of equations for determining 

We shall suppose that the coefficients a are of such sort that each of the 
expanded determinants in (3) has explicitly present at least one of the 
parameters A. Then system (3) is a set of algebraic equations r in number 
and involving r unknown quantities. We assume furthermore that the 
coefficients a are such that these r relations are independent in such wise 
that we have only a finite number of sets of characteristic values. 

These sets of characteristic values we shall then denote by 


for varying values of p, the two ordered sets being distinct for two distinct 
values of p. The corresponding solutions of (1) and (2) we shall then denote 
by 
3. Conjugate Character of the Solutions of (1) and (2). With the notation 
adopted in the preceding section for solutions of (1) and (2) we have relations 


which may be written in the form 


Nh 


(5) (Aonij -+ + NP arnij) xi?) = 0, aj 2, 
j=l 
Nh 
(6) Dy (Gongs + + 
j=1 
nm, 
— = 0, t= 1,2, 
j=i e=l 


systems being formed for each value h of the set 1, 2, ---, r. For fixed h 
and 2 let us multiply the former of these equations member by member 
by y‘? and the latter by — 2‘), and let us then add the two resulting equa- 
tions member by member. Still holding h fixed in the relation so obtained, 
sum as to 7 from 1 to n,. That part of the first member of the resulting 
equation which comes from the first member of (5) and the first line of the 
first member of (6) balances to zero and we have a relation which may 
readily be reduced to the form 


r nh Nh 


(7) — r”) = 0, 1, 2, «++, 
= j= 


s=1 
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Here we have r equations involving the r quantities NM” — \, s = 1, 2, 
, r, these being not simultaneously zero when p and o are different. 
Hence the determinant of the coefficients of these quantities in (7) must 
have the value zero. If we employ different summation variables 7 and j 
for the different rows in this determinant, the relation may be put in the 


form 


| nN», Nn 


p) 
= 0, p + 
| 


where the element written explicitly is that in the sth column and the hth 
row. It is easy to see that this relation may be put into the following more 


convenient form: 


where D,,;,...;,;, denotes the determinant 
| 74; 213, %, i 


_ | °° * 


| Qirj,i, Qerj i * Urrj,i, | 


This relation for varying p and oa, expresses the fundamental conjugate 
character of the solutions of (1) and (2) with respect to each other. 

In (3) let the Ai, ---, A, be replaced by the set of characteristic values 
AP, --+, A”. For a given h let a’, denote the cofactor of the element in 
the 7th row and jth column of the determinant in the resulting first member. 
Then, if for each h some one of these cofactors is different from zero, we 


have a relation of the form 


Similarly, we have 


Employing &;, and 7; as the (necessarily non-zero) factors of proportionality, 


we have 
— (p) (p) — (p) 


If we denote the first member of (8) by S(p, «) we have then 


| 

D,;, 
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Now the second member of (4) depends only on the coefficients in equa- 
tions (1) and (2). If 7 and j exist so that this second member is different 
from zero we have S(p, p) + 0. 

Since we shall be interested in these algebraic results primarily for their 
use as a heuristic guide in certain transcendental problems we shall naturally 
exclude from consideration such exceptional cases as are complicated or 
otherwise fail to serve our purpose. Consequently we shall suppose now 
and usually henceforth that the coefficients of (1) and (2) are of such character 
that the second member of (9) is different from zero for every p. Then we have 
the following result: 

The first member of (8) vs equal to zero when p + o and is different from 
zero when p = oa. 

In the special case in which (1) is self-adjoint and the 2’s and the y’s, 
for a given p, are taken to be the same, it is clear that the condition that all 
determinants in (8) shall have the same sign (excluding therefore zero values) 
is a sufficient condition for the non-vanishing of the first member of (8) 
when p = oa. 

In order to bring out clearly the nature of the conditions of conjugacy 
which we have just established let us consider the special case in which 
r= 1. Systems (1) and (2) may then be conveniently written in the form 


(10) (ci = 0, = + = 0, 2, 
j=1 j=l 
The conditions of conjugacy reduce to 


= Oife + p 
(a) 
(11) bit Y; + 0ifo= 
If we have the more special case in which );; = b;, b;; = 0 if J + 1, these 
conditions become 


(p),(o)) Oif o + p, 
(12) bia Py 
When b; = 1 for every 2 these become merely the usual conditions of bi- 
orthogonality; and these in turn reduce to the usual conditions of orthogon- 
ality in case the two systems in (10) are identical and the same solutions 
x and y are taken for the two systems. From this it is seen that our formule 
give extensive generalizations of classic elementary relations of wide use- 
fulness. 

Let us consider a special case of equations (1) and (2) giving a rather 
simple generalization of the results in the preceding paragraph. Let us 
suppose that each system in (1) contains but a single parameter, say that Aj, 
is the only parameter in the system for a given value of h. Then the r 


| 
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equations (3) involve each a single one of the r parameters \i, ---, A,, so 
that the possible values of any given parameter are obtained by solving a 
single equation explicitly given. All possible combinations of solutions 
are allowable in forming the various sets of characteristic values. The 
determinants in (8) reduce in the present case to their diagonal elements, 
so that we have for this case the simpler conditions of conjugacy 


»\{ =0ito +p, 
(13) ) | + Oifc = p. 


4=1 i,=1 
If for a given subscript / no characteristic value of \; is equal to zero 
we may divide systems (1) and (2) by \; and so introduce a new set of r 
parameters 1/A,, A,/A, for i + k. These enter linearly and in the same 
way as the original parameters. It is easy to see that the solutions x and 
y are the same as before. Consequently we have new sets of relations 
similar to (8) and (13) and differing from them only in having a,,;; re- 
placed throughout by 
4, Expansions by Use of Conjugate Sets.—Let m be the number of sets 
of characteristic values of \i, ---, A, and let us suppose that the constants 
; for 2, varying from 1 to mn for each h may be “expanded”’ in the form 


c=] 

where the c; are independent of the subscripts 7;.. Then the properties of 
conjugacy developed in § 3 are available for an immediate determination 
of the values of the c,—much after the manner of determining the coeffi- 
cients of a Fourier ae For this purpose we multiply both sides of 
equation (14) by y{),--- times the determinant D,,;,...;;, explicitly 
written in connection w ith (8) and in the resulting equation sum as to the 
i and the 7. In view of the conditions of conjugacy the resulting quantity 
in the second member depends on only one of the constants ¢c;,, namely ¢;; 


and we then have 


ny n Np Ry 
1=1 =1 j,=1 


ji 


It is easy to obtain broad sufficient conditions for the validity of the 
“‘expansion”’ (14) for any given set 2;,;,...:, of mim2 +--+ m, constants. We 
suppose that m = nyn2 --- n, So that the number of distinct solutions of (3) 
is equal to the product of the degrees of the several equations in that system 
If then we look upon (14) as a system of m equations for determining 


{ 
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the quantities ¢1, C2, +++, Cm it is sufficient that the determinant D of this 
system shall be different from zero. But we have seen (in §3) that 
a) = £,a, for some appropriate 7 (dependent perhaps upon k), the 
quantities £, being different from zero. Hence the determinant D will be 
different from zero if the corresponding determinant formed with a‘, in 
place of 2%? be different from zero. This last is a condition on the coeff- 
cients a in equation (1). Thus we may express in terms of the coefficients 
in (1) broad general conditions which are sufficient to ensure the validity 
of the expansion (14), where the coefficients c; have the values given in (15). 

It is not in accordance with our present plan to seek those modifications 
of the foregoing expansion which are necessary in the exceptional cases 
excluded from the preceding discussion. It does not appear that we shall 
need the theory for these cases in order to avail ourselves of the heuristic 
guide afforded by the algebraic theory in its non-exceptional and more 
general aspects. 

In order to bring out clearly the essential simplicity of these formule 
let us observe that for the special cases involved in equations (11) and (12) 


we have 


n 


k=] 


where ¢; has one or the other of the values 


n n 


n 


j= 


= n n 


i=1 j=l 


according as we have the case of (11) or of (12). Ifb; = 1 the latter reaches 
the maximum of elegance, at least if equations (10) are then self-adjoint 
and the same solutions x and y of the two systems are taken. 

For the general case when m = nin2 --- n, there is a certain gain in 
suggestiveness if we denote the parameters and solutions of (1) and (2) by 


(17) Ck 


(ky ke... k,) (ky 


where k, runs over the set 1, 2, ---, ms. Then (14) takes the form 


ny 


(18) = >» TT »), 


The coefficients c are determined in the same manner as before. Just as we 
may look upon (16) and (17) as the analogue of expansions of functions of 
one variable in orthogonal functions of that variable (including the special 
case of Fourier series) so shall we be able to look upon (18) as the analogue 


| 
| 
| 
ij 


78 CARMICHAEL: Boundary Value and Expansion Problems. 


of like expansions (including double Fourier series) by which functions of 
several variables are expressed in terms of products of functions each of one 
variable—such expansions being of the kind used by the phsyicists in the 
solution of many classical problems. Indeed these are much more than 
analogues; transcendental expansion problems of several well-known types 
(including those mentioned) and of many new types are as a matter of fact 
limiting cases of the algebraic expansion problem which we are treating 
(and we are developing the latter primarily for its applications to the former). 

In the foregoing expansion formule the number of subscripts on the 
“function” z to be expanded in terms of a given set of constants by aid 
of the property of conjugacy is equal to the number of parameters A involved 
in the original system of algebraic equations. When we proceed to our 
transcendental limiting cases we shall see that these algebraic results are 
in the form best suited to applications to equations involving functions of 
one variable, as ordinary differential equations, ordinary difference equa- 
tions, integral equations, and certain types of mixed systems. In order to 
obtain a form suitable as the heuristic guide in problems involving partial 
differential equations or equations of any type in which we have functions 
of more than one variable in the original equations, it is desirable to look 
upon our expansion formule in a way slightly different from that in evidence 
in the foregoing work. 

We can best bring out what is needful by considering first the particular 
expansion of equation (16). Let us suppose that n is the product uv of two 
integers. Let us replace the subscript 2, running over the set 1, 2, ---, n, 
by the double subscript 7j where 2 runs over 1, 2, ---, wu, and j over 1, 2, 
--+,yv. Then equations (16) become 


The coefficients c; are of course determined in the same way as before but 
are naturally expressed in a form differing from that in (17) by having the 
single summation as to 2 replaced by a double summation as to 2 and j. 

5. Graphic Representations of Solutions.—Let us consider the graphical 
representation of the set 1, w2, «++, Un of constants in the following manner. 
On any given straight line segment (as for instance the points s such that 
a = s Sb) let us erect n perpendiculars two of which are at the ends of the 
segment while the other n — 2 are evenly or evenly distributed on the 
interior of the interval. Let these be marked from left to right by the 
numbers 1, 2, ---, m; and consider them as analogous to the n codrdinate 
axes of a space of n dimensions. Considering (21, we, +++, Un) aS a point in 
space of n dimensions we may now represent this single point by a certain 
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junction of n points in our figure, one on each of the axes, the 7th one being 
on the 7th axis and at a distance |w;| from the original segment and above 
it or below it according as wu; is positive or negative. Now join by straight 
line segments the point on each interior axis to the points on the two adja- 
cent axes. We thus obtain a broken line. We shall say that this broken 
line is the graphic representation* of the point (a2 --- %,) in n dimen- 
sions or of the set of constants w, w2, +++, Un. The points at which (values 
of s for which) this broken line cuts the original line segment (viewed as the 
axis of s) we shall call the zeros of the set w;, in analogy with the usual ter- 
minology for the zeros of a function w/(s). 

Now let us suppose that the same set of axes is used for the representa- 
tion of the n sets x, k = 1, 2, ---, n, and the set z; of equation (16). 
Let x (s), k = 1, 2, ---, n, and 2(s) denote the functions of the continuous 
variable s represented by the broken lines corresponding as indicated above 
to the foregoing sets of constants. Then we shall show that 


where the constants c; have the same values as in (16). For this purpose 
it is sufficient to consider the graph between two consecutive vertical 
axes, say between the ith and the (c+ 1)th. Let s; and s;,; be the values 
of s at the intersections of these vertical axes with the s-axis. Then between 
the two vertical axes in consideration we have 
Eom 1, 2, 
From these relations and from (16) with the subscripts 7 and 7 + 1 it follows 
that (20) is valid for the interval in consideration. Hence (20) is valid 
for all values of s on the range a = s = b of the original line segment. 

Let us consider the like matter for functions of two subscripts and 
expansions of the form (19). For the representation of w;; for 7 = 1, 2, 
-++,u,andj = 1, 2, ---,», we shall start from v planes, one for each value of 
j. In each of these we use for s the same range a = s S b and place the 
planes in order 1, 2, ---, v for j each directly in front of the preceding one 
and arrange the vertical axes so that the axes in any one plane are the ortho- 
gonal projections on that plane of the axes in any other plane. Then for 
fixed 7 we extend the w,;; to wu; (s) by linear interpolation as before. Then 
for each value of s we connect the points w;(s), 7 = 1, 2, ---, v, by straight 


* Essentially the same graphic representation has been employed by M. B. Porter in 
Annals of Mathematics (2), 3 (1901), p. 56. 
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line segments joining the consecutive points. We thus get a sort of broken 
surface affording the graphical representation of a function u(s, ¢), gotten 
(we may say) by linear interpolation from the points w;;._ We call this the 
graphical representation of the set of constants u;; for i = 1, 2, ---, wu, and 
j= 1,2,---,v. Lete St =d define the range of ¢ in this representation, 
this range evidently depending on the positions of the v planes employed 
in setting up the graphical representation. 

Let us replace the subscripts 7, 7 in the functions in (19) by the con- 
tinuous variables s, ¢ in accordance with the method just indicated. It is 
easy to prove that we then have 


(21) 2(s,t) = Doew(s,t), asssb, 
k=1 


where the constants c; have the same values as in equation (19). In fact 
the argument associated with (20) shows that we have 


2j(s) = j=1,2,---,». 


A repetition of the same argument for these last equations with s fixed 
leads at once to relation (21). 

Similar extensions may be associated with equation (14) so that we 
readily come through to the relation 


m 
2(81, Sr) = (81) 22 (S2) 
k=1 


where the c; have the same values as in (14) and the functions involved 
depend on continuous variables. Further generalizations may also be 
made in the direction suggested by the last paragraph of § 4. 

6. On the Distribution of Zeros of Solutions.—The graphic representation 
of solutions described in the preceding section we shall find to afford one 
of the most useful ideas in the heuristic guide to transcendental problems 
arising from the algebraic theory which we are now developing. In order 
to have at hand an aspect of its usefulness not brought out directly by the 
applications to expansion problems, let us consider a certain fundamental 
question concerning the relative distribution of zeros of the solutions of 
certain types of algebraic systems. 

Let us first consider a system of n homogeneous linear algebraic equations 


n+2 
(22) i=1,2,---,n, 
j=l 
in the n + 2 unknown quantities 2, 22, +++, %n42; and let us suppose that 


the matrix of this system is of rank n. The system has two solutions x; = u; 


| 
| 

| 
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and x; = v; such that every solution is expressible in the form x; = cu; + dv,;, 
where ¢ and d are constants, while conversely every such 2; is a solution. 
The sets wu; and »; are linearly independent. If 2; = u;and 2; = 0; are any 
two linearly independent solutions and if A;,; denotes the determinant of 
the matrix obtained from the matrix of coefficients in (22) by striking out 
the kth and /th columns, it may be shown that Aj; and the determinant 
Uzv1 — Uw; are both zero or neither zero, provided that when A;, = 0 we do not 
have Am: = 0 = Axm for every m of the set 1, 2, ---, n — 2, except k and l. 

In order to prove the last statement, let us write equations (22) in the 


form 


n+2 
j= 
Then we have 
(23) Akitm =+ k A km%1, 


the ambiguous signs + being obviously determinate. 

Let us suppose that A;; + 0. Then we can not have simultaneously 
uy; = 0, u, = 0; for, if so, (23) would imply that w; = 0 for every2. Like- 
wise we can not have simultaneously v; = 0, v7, = 0. Let wu, be one of the 
two quantities u,, uw: which is different from zero. Then if we suppose 
that u,v; — uw; = 0 we have vy = ru, and vo; = ru; for an appropriate 
constant r; whence it would follow through use of (23) that v; = ru; for 
every 2, so that the solutions wu; and v; would be dependent, contrary to 
hypothesis. Therefore, if A;: + 0 we have u,v; — uw, + 0. 

Suppose next that Ay; = 0. Substituting in (23) the solutions wu; and 
v;, we have 

+ Anite + Akmui = 0, + + = 0. 
Then, if wxv, — uw, + 0, we have An: = 0 = Agm for every m of the set 
1, 2, ---, n+ 2, except k and /, contrary to hypothesis. Hence, when 
Axi = 0 we have — Uv, = O. 

This completes the proof of the statement in consideration. 

The foregoing result may readily be extended to the case of a system of 
rank n of n equations in n + h unknown quantities, where h is any integer 
not less than 2. If a linearly independent set of solutions is w1;, U2;, + -, 
ux; We compare the determinant |wjk;| whose element in jth row and 7th 
column is u;k; with the determinant Aj,;,...;, of the matrix of order n 
remaining when one strikes from the matrix of coefficients of the original 
system the columns with numbers fi, ky, ---, ky. For a given set ki, ke, 
-++, k, these determinants are both zero or neither zero, under conditions 
similar to those named in the preceding case. 

If we replace x in (23) by u and by v we obtain readily from the resulting 


q 

| 

a 
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equations the relations 


Um Ul UR Ul 
(24) Axi =+ Ami 
VI UK VI 
and, in particular, the relations 
| 
| Uk+1 | UK 
| Uk+2 | Vk Uk+1 


Hence, if A;, +; for a given range of consecutive values of the integer is 7 
of one sign (excluding the value zero), then the determinant w,;, 


is of one sign for every 2 of the same range. It is obvious that this result 
also may be appropriately extended to the more general case mentioned in 
the preceding paragraph. 

Let us now prove the following theorem: 

Suppose that A,, ;4: for a given range R of consecutive values of the 
integer 2 is of one sign (and hence not zero). Let k and / be two integers of 
this range, and suppose that uw, and w;41 are not of the same sign and that 
if either of them is zero it is wx (in which case u;4; is certainly different from 
zero). Suppose again that w; and w:,; are not of the same sign, and that 
/ is the smallest value of the subscript greater than k for which this is so. 
Then an integer h exists, k = h =1, such that v and v4; are not of the 
same sign. Moreover, there are not two such values of h except possibly 
for the pair h = k, h = 1. 

The last statement follows readily from the preceding part of the 
paragraph; for, if it were not true, the earlier part of the theorem could 
be applied with uw and »v interchanged so as to show that the consecutive 
character of u; and uw; is not maintained. 

For the given range R of values of 7 the determinant w; is of one sign, 
as we have already seen. Without loss of generality we may suppose that 
w; > 0, since if it were less than zero we should merely have to change the 
sign of every u; or of every v7; Then we have 

> Uj410j. 

Let us first consider the case in which u; = 0. Then v; is of the opposite 
sign to uz41. Without loss of generality in argument we may (and we will) 
take positive and v; negative. Then ---, ware all positive while 
W141 1s zero or negative. If v, is zero or negative we see from the inequality 
> Ur410, that is positive. Hence either or is positive. 
Hence there is a change of sign in the sequence 2x, 441, -+ +, 0241 SO that the 
statement in consideration is true for the case when u; = 0. 


| 
| | 
: wi = | ly 
Vv; Vi+1 
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If uw, + O we may take it to be negative. Then wz41, ---, uw, are positive 
and is zero or negative. From the inequality > We see 
that if v; is positive then v,41 is negative and the truth of the statement in 
consideration is granted. It remains to consider the case when 2; is nega- 
tive. Either the statement is verified or we have vj41, ---, v; all negative. 
In the latter case we see from the relation wyvi4; > ui41v; that the second 
member is zero or positive and hence that v1 is positive. Hence v; and 
0141 are of different signs and the statement is verified in this case. 

This completes the proof of the statement in consideration. 

Let us now consider the relative distribution of the zeros of the functions 
u(s) and v(s) gotten from the set of constants wu; and v; by linear interpolation 
after the method of §5. As in the immediately preceding discussion we 
confine attention to a sequence of consecutive intervals corresponding to a 
range R of consecutive values of the subscript 7 for which Aj, ;:; is of one 
sign. From the results just proved in the foregoing paragraphs it follows 
that there is a zero of v(s) between two consecutive zeros of u(s) except 
possibly for the case when the h of these paragraphs has the value k or /. 

Let us consider the case when h = k. We have 


= UK+1Vky 


while v,41 and v; are not of the same sign and are not both zero. If u; = 0 
we have a case when v(s) has a zero between the two given consecutive 
zeros of u(s). Then consider the case when u; + 0. Then wz and wz4 
are of different signs; and we may take the former to be negative and the 
latter positive without loss of generality in argument, and this we do. If 
in the interval joining the points corresponding to k and k + 1 the zero of 
v(s) is to the right of that of w(s) we have a zero of v(s) between the two 
consecutive zeros of u(s). We then consider further the other case, namely, 
that in which the zero of v(s) in the interval in question is to the left of 
that of u(s). In this case v4.41 + 0. If v4.41 is positive we have 


ke 


which is contrary to the present hypothesis that the zero of v(s) in the inter- 
val in question is to the left of that of u(s). Hence we must now have 
Vegi negative. Then either v(s) vanishes between the two: consecutive 
zeros of u(s) or v7; is negative. We consider the latter possibility. Now w, 
is positive. Hence from the relation wi41 > wi4101 we have 


v1 Ul 
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From this and the fact that wu; and w;,; are not of the same sign it follows 
readily that the zero of v(s) in the interval from / to 1 + 1 is to the left of 
that of u(s) in the interval, and hence in this case that a zero of v(s) occurs 
between two consecutive zeros of u(s). We have now exhausted all the 
possibilities for the case when h = k and have found that in this case a 
zero of v(s) certainly occurs between two consecutive zeros of u(s). 

The case when h = / may be treated in a similar way, with the same 
conclusion. 

From this it follows that there is always at least one zero of v(s) between 
two consecutive zeros of u(s) on the interval in consideration. There can 
not be more than one; for, if so, the result could be applied with w and » 
interchanged so as to lead to the conclusion that u(s) should have a zero 
between two consecutive zeros of the same function. 

The result thus obtained may be formulated into the following theorem: 

THEOREM. Let Aj, i41 for a given range R of consecutive values of the 
integer 1 be of one sign and let I denote the interval of the s-axis corresponding 
to this range of 1 in the sense of the treatment in § 5. Let u; and v; be two linearly 
independent solutions of system (22); and let these solutions be extended by the 
method of linear interpolation employed in § 5 to the functions u(s) and v(s). 
Then on the interval I the zeros of u(s) and v(s) separate each other. 

Let us now consider more generally a system of equations 


n+h 
in the n+ h unknown quantities 7, %2, +++, in, where h = 2; and let 


us suppose that the matrix of this system is of rank n. The system then 
has h linearly independent solutions; we designate such a set by 


h 
xs (h) 


) 


Let D; be the determinant of order n of the square matrix gotten from 
the matrix of coefficients in (25) by striking out h consecutive columns 
beginning with the 7th; let R be a range of consecutive values of the subscript 
i for which D; is of one sign (zero values being excluded); and let J denote 
the interval of the s-axis corresponding to this range of 7 in the sense of 


the treatment in § 5. 
Now let w; denote the determinant 


| | 
2) (2) | h—1,,(2) 
h Ah) Ah | (h Ah —1,,(/) | 


where A, A’, ---, A’, in the last form of the determinant, denote the first, 


— 
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second, ---, (h — 1)th differences of the x; with respect to the variable 7. 
It is obvious that the two forms of the determinant have the same value. 

Since D; is of one sign for 2 on R it follows from the earlier work in this 
section that w; is of one sign fori on R. Let us consider the matrix of any 
m rows of the first [second] determinant form for w; in (26) and let 
Aimjl_Cimj |, j = 1, 2, +++, jm, be all the m-rowed determinants formed from 
this matrix. Let Bim;_Dimj],j = 1, 2, +--+, jm, be the corresponding alge- 
braic complements. By means of the Laplace development for determi- 
nants we then have 

Jm Jm 
(27) ws = AimjBimj = CimjD 
j=l j=l 

Suppose that k is a value of 7 on the range R such that Bimj [Dimj_] is 
of one sign or zero for each j of the set 7 = 1, 2, ---, jm; and suppose that 
lis a value of 2 also on R for which these quantites are of the opposite sign 
or zero. Then, since w; is of one sign for 2 on R, it can not be true that all 
of the quantities A imj, Aim; [Ckmj, Cimj_| are of one sign for the fixed value 
of m and for j varying over its whole range or indeed over just that part 
of its range for which the corresponding Bim; or Bimj Dimj or Dimj| has a 
value different from zero. 

In particular, if Bim; and and Dim; both have the value 
zero except for a single value p of j and for this value have different signs, 
then Aim, and Aimy [Cim, and Cim,] have opposite signs. Hence the 
function Asmp [Csm,_] obtained from A im, ['Cimp_] by linear interpolation as in 
§ 5 has a zero on that part of the s-axis which corresponds to the range of 2 
from k to / inclusive. 

Let us now consider the first determinant in (26) for the case when m is 
unity and the one row singled out is the last one. Then the Aji;, for 
varying j, are the elements of the last row of the determinant and the Bj; 
are their cofactors. If for 1 = k these cofactors are of one sign or zero and 
for 7 = 1 are of the opposite sign or zero, where k and / are both on the 
range R of 7, then not all of the quantities 2%, af}1, --+, zf2,_1 are of one 
sign. Hence if we form the function x (s) by linear interpolation from 
a” this function z(s) has at least one zero in the interval corresponding 
to the range of i from k tol-++ h—1. Thus for the cases when our results 
apply we are able by means of h — 1 linearly independent solutions of (25) 
to define an interval for 7 and a corresponding interval on the s-axis such 
that every solution 2z{ of (25) which is linearly independent of the first 
h — 1 solutions gives rise by linear interpolation to a function «™(s) which 
vanishes on the interval of the s-axis in question. It will be observed that 
this result is a direct generalization of a portion of that obtained above for 
the case when h = 2. 
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Let us consider the second determinant in (26) for the case when m 
is unity and the one row singled out is the last one. Then the (;, for 
varying j, are the elements of the last row and the D,; are their cofactors. 
If the latter for 2 = k are of one sign or zero and for i = | are of the opposite 
sign or zero, where k and / are values of 7 on the range R, then the quantities 
a, Ax®, «++, Aba”) for the two values k and / of 7, can not all be of the 
same sign. In particular, if they are all positive for 7 = k, then one at least 
of them must fail to be positive at 2 = /; so that at least one of them changes 
sign or becomes zero in the interval in question. 

We may look upon the modification of the first determinant in (26) to 
produce the second as a single instance of a great variety of transformations 
of that first determinant, each being made by replacing columns in that 
determinant by various linear homogeneous combinations of columns (the 
coefficients being constants or functions of 7) in such wise that the non- 
vanishing of w; on a range R implies the non-vanishing of the modified deter- 
minant w; on the same range R. For every such w; we may obtain a se- 
quence of theorems of the sort just brought to light in the treatment of the 
two forms given in (26). Such properties as are thus indicated will be 
useful to us later in the study of certain transcendental problems which are 
limiting cases of the algebraic problems now in consideration. 

Solutions of (25), for a given value of h > 2, may also be studied by 
means of similar systems for a smaller value of h, this smaller value of h 
being not less than 2. For this purpose we adjoin to (25) one or more 
additional equations in the same unknown quantities, these equations being 
linearly independent of each other and of those in (25). In particular, the 
system may thus be replaced by a new system of the same form but with 
only two linearly independent solutions. The foregoing results for such a 
system, when applied to this, yield properties of those solutions of (25) 
which satisfy the auxiliary conditions, the latter being somewhat analogous 
to boundary conditions for differential equations. A considerable hold 
can thus be gotten on a class of solutions of (25). This obvious remark will 
serve a useful heuristic purpose in studying certain properties of differential 
equations with boundary conditions. 

In this section we have confined ourselves principally to that part of the 
range of 7 for which D; is of one sign. We propose the problem of investigat- 
ing the matter for the whole range of 7 from 1 to n and thus of associating the 
zeros of D; and of certain combinations of the solutions of (25) for the whole 
range of 2 from 1 to n. - 

7. General Character of a Fundamental Limiting Operation.—It is well 
known, as we have already pointed out more fully, that several important 
types of functional equations, such as linear differential equations, linear 
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integral equations, etc., can be exhibited as limiting cases of algebraic 
systems. In a large range of limiting cases the conjugate character of the 
solutions of (1) and (2) which we have discussed in § 3 and the expansion by 
use of conjugate sets treated in § 4 are carried over to the limiting problems 
in such a way as to afford a definite suggestion of properties to be expected 
and of methods of analysis for examining them. We shall now consider 
the general character of the limiting operation in question with reference 
to the properties just indicated. 

The graphic representations employed in § 5 afford a starting-point. 
Let us take a fixed interval a = s = b of the real s-axis and let us associate 
with certain points of this interval the discrete values of 7 and 7 in equations 
(1) and (2). For the Ath equation in either system we take on (ab) a set of 
points n, in number including the points a and b; and we interpolate x; 
into a function 2,(s) by linear interpolation as in $5. Similarly, a function 
Ani; iS interpolated into a function a;,(s, t) by a method also described in § 5. 
“quations (1) and (2) may now be looked upon as establishing relations 
among the functional values of these functions of s and ¢ at certain points 
only of the axes of s and t; and so of defining the solution functions at these 
points and these alone, their definitions being completed by the method of 
interpolation agreed upon. 

The limiting processes described consist essentially in this: let the various 
numbers n,, (or at least a part of them) increase indefinitely; and as they 
increase let the functions 2,(s) pass through a corresponding sequence of 
changes. Various limiting processes are thus set up in the original equations 
and in their solutions. If these lead to replacing the original equations by 
well-defined equations and their solutions by well-defined functions, we 
have in the process a suggestion of a heuristic guide to probable solutions 
of the limiting problems and to certain probable fundamental properties 
of these. As we proceed it will be seen that this guide has great usefulness 
in a considerable range of problems. 

We shall usually require that the distribution of basic points on the 
interval (ab) of the s-axis shall undergo change in such wise that the maxi- 
mum distance between two consecutive points shall approach zero as np, 
becomes infinite. This suggests that the sum as to 7 in (1) and (2) may 
sometimes be replaced by an integral as to ¢ from a to b; and, when this is 
to be done, it is desirable that the coefficients a;,,;; shall be replaced each by 
a product one factor of which is to correspond to the differential in the 
limiting case. [If the basis points on the s-axis (or f-axis) are equally 
spaced one may introduce this needed factor simply by multiplying the 
equation through by As (or At), the distance between two consecutive basic 
points. This is perhaps the easiest way to make the modification in ques- 
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tion. ] In fact such a form has frequently been employed in setting up a 
limiting process to guide one to the theory of integral equations, first by 
Volterra and later by others. In numerous ways (of which that just men- 
tioned is one of the most important) one may subject the coefficients in 
(1) and (2) to subsidiary conditions such that the limiting process in con- 
sideration shall lead to desired forms of functional equations as limiting 
cases, whether differential equations or integral equations or equations of 
other type. 

The several systems in a single set of systems (1) or (2) are independent 
of each other in such wise that the limiting processes may be applied to 
each with a certain measure of independence of the others, so that we might 
have for instance an ordinary differential equation and an integral equation 
to make up the limiting system of two equations for the case of two param- 
eters \. Numerous other combinations are logically possible; and several 
of them are interesting for the results to which they give rise. In a second 
paper we shall give a precise characterization of several problems to which 
we are led in carrying out the limiting operation in a variety of ways or 
combinations of ways. 

It remains here to note certain characteristic results for the finite case 
which persist (in modified form of course) in a great variety of limiting cases. 

The property of conjugacy expressed by the equations of system (8) is 
an important one which persists after certain of the limiting processes just 
indicated have been carried out. It takes a large number of forms owing 
to the great variety of ways in which the limiting process may be set up. 
One of the simplest and at the same time most useful of these is that in 
- which the determinant D;,,;,...:,;, in (8) becomes a function D(s:, 
s,, t,) and the summation in (8) is replaced by an integration, the algebraic 
system having first been taken in a form which is convenient for passing 
to the limit in a manner to realize such a result. The property of conjugacy 
may then be expressed by the relations 


b b b r 
(28) f D(s1, th, +++, Sr, tr) II (sn) -ds,dt;- - -ds,dt,=0, 
a va a 
p +o. 


For a wide range of cases we have in the limiting problem an infinite 
number of sets of characteristic values for the parameters \ and hence an 
infinite number of sets of functions xs), y(t), h = 1, ---, r, for varying 
p, for which relation (28) is valid. 

One of the most frequently occurring and important special cases of 


(28) is that which corresponds to equation (12), namely, 
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When we proceed in the second paper to certain important instances 
of the theory here sketched roughly in a general sort of way we shall find 
that equations essentially like (28) or special cases of it appear in a great 
variety of forms; but a little examination of each system will show that 
its equations arise in a natural way as limiting cases of the equations in 
(8) and that one is led to them by the guide afforded by the general con- 
siderations barely outlined here. 

Now if the repeated integral in the first member of (28) has a value 
different from zero when p = o (a condition which is realized in a wide 
range of important cases, as we shall see) and if a function f(s1, s2, ---, $,) 
has an expansion of the form 


f (81, 82, 8) = 


where the ¢, are constants and the values of p correspond to the various 
sets of characteristic values of the parameters \, then the coefficients ¢, 
are readily found by a method analogous to that by which the Fourier 
coefficients for a function are found; and we have in fact 


b b b 
i) (81, 82, +++, 8)D(s1, th, +++, Sr, te (te) 


-- - ds,dt, 


b b b 
D(81, th, Sry &) II ap) t,)-ds,dt, --+ ds,dt, 


h=1 


In each instance of our theory this formal expansion or the requisite modi- 
fication of it is used to set the basic expansion problem associated with the 
instance in consideration. 

In the foregoing discussion the number of variables in the function f to 
be expanded is equal to the number of parameters in the algebraic problem 
from which we passed to the transcendental problem. This correspondence 
is not essential to the nature of the process involved. The desired extension 
‘an best be brought out by starting from the particular expansion (19). 
We may carry out our method of linear interpolation as at the close of § 5 
and obtain the expansion (21). If now we proceed to the limiting case in 
such wise that » and v simultaneously approach infinity and if we require 
that the distribution of axis points on the s-axis and the f-axis shall have the 
usual property of condensation and if the number of distinct characteristic 
values approaches infinity during the process, we shall be led heuristically 
to an expansion of the form 


‘(s, t) = ¥ cpr (8, t), 


k=1 


= 
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where the 2’s are now solutions of the limiting problem. It is easy to see 
that the properties of conjugacy are maintained (under appropriate cireum- 
stances) and that formule are therefore easily found for the values of the c,. 
We do. not attempt to investigate in general the range of validity of these 
results. We shall later find in them a useful heuristic guide to problems in 
partial differential equations; and for the cases involved in such applica- 
tions we shall analyze the matter of range of validity. 

In a similar way we may obtain numerous modifications of the foregoing 
expansion of f(s1, so, ---, 8.) by starting from the modifications of equations 
(14) suggested by the method of the preceding paragraph for the simplest 
special case and proceeding in a similar manner. 

Another method of proceeding from (19) to a limiting case will be very 
useful to us in dealing with a system of ordinary linear differential equations. 
In this method we suppose that the passage to the limit is effected in such 
way that w remains fixed so that 7 is always a discrete variable ranging 
over the finite set 1, 2, ---, u, while v becomes infinite and the variable j 
passes in our usual way to the continuous variable s in the limiting equation. 
Then expansion (19) passes over into the form 


fils) = Does), 

a form which may be looked upon as affording the simultaneous expansion 
of w given functions f;(s) of a continuous variable s in terms of u sets of 
functions 2;,(s), the coefficients c; being the same in each of the two expan- 
sions.* Obviously, numerous more general forms of like expansions of 
systems of functions of one or more variables may be obtained by passing 
from the general relation (14) to limiting cases in the various ways sufh- 
ciently suggested in the foregoing instances. 

The details of the matters suggested in the three preceding paragraphs 
will be developed only in connection with the instances of certain tran- 
scendental problems to which they afford the requisite guide. 

It is doubtless true that many other properties of the algebraic systems 
can be carried over to certain of the important limiting cases. In fact 
certain of the zero separation theorems of § 6 can be utilized to suggest like 
theorems for the solutions of certain transcendental problems. Moreover, 
the interaction of the two types of problems may well be made to work 
in both directions. A part of the results in § 6 were indeed suggested by 
classic properties of the solutions of linear differential equations of the 


* A special case of expansions of this sort first came to my attention in the dissertation 
of Dr. C. C. Camp (a manuscript copy of which he was kind enough to lend me); and the 
suggestion afforded by it has been very useful to me in the present investigation. 
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second order. Once in hand for the simpler cases these results were easily 
extended to other cases of algebraic systems and then afforded a new guide 
to certain related transcendental problems. There is thus a repeated 
mutual reaction between the two types of problems which is useful in the 
development of the theory of each. 

8. A Second Type of Limiting Case.—Let us now consider the r homo- 
geneous linear systems in the infinite number of unknown quantities 2, 


(29) (Qonij + -+- + = 0, = a, 


a separate system being formed for each value h of the set 1, 2, ---, 1. 
With these systems let us associate the r adjoint systems 


8 
(50) (Gongs + + = O, t= 1,2, 


Here \i, Ae, «++, A, are 7 parameters by the choice of which the simul- 
taneous consistency of the r systems in either set is to be secured. For these 
infinite systems a great variety of possibilities may be realized by special 
methods of defining the coefficients a; and a comprehensive analysis of all 
mutual possibilities would require the examination of a great many cases. 
Consequently we shall confine our attention to a range of cases in which the 
theory may be set forth in the simplest manner. 

A set of values of \u, ---, A, for which (29) [((30) ] has a solution’ not 
identically zero will be called a set of characteristic values for (29) [(30) ]. 
In the case of the finite problem of § 2 we saw that the sets of characteristic 
values for the two adjoint problems were identical. For the infinite case 
of this section we shall consider only those systems for which this property 
is maintained. It is clearly maintained in cases in which those character- 
istic values can be determined by determinantal equations analdgous to 
equations (3) of § 2. If this method of determining the sets of characteristic 
values is not suitable we may proceed as follows for a certain range of cases 
of equations (29) and (30); omit from each system r equations, one for each 
value of h; solve the remaining equations as non-homogeneous equations, 
after having assigned a value to one of the unknown quantities (an operation 
which can be carried out by known processes in a wide range of cases), and 
suppose that this solution is unique (as indeed it is in a wide range of known 
cases); now substitute into the equation omitted from each system the 
values of the unknown quantities so obtained; thus we have for system (29), 
and also for system (30), r equations involving the r parameters; their 
solutions give the sets of characteristic values. 

Now, we shall suppose that the sets of characteristic values, however 
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determined (whether by one of the methods indicated or by another), is 
denumerably infinite and we shall denote them and the corresponding 
solutions by 
20, p= 1, 2,8, 

Several of the properties developed for the corresponding algebraic problem 
are now readily carried over to these solutions of a limiting problem, at 
least in a formal way. Since our present interest in the general procedure 
arises from intended applications of it to relatively simple instances we 
shall not undertake an analysis of the range of validity of our results, but 
shall content ourselves with little more than an exhibition of the formal 
processes. 

If we proceed by a method in all respects similar to that employed in 
§ 3 and if we have before us a case in which the order of infinite summations 
may be interchanged in ways now needful in carrying out the same method, 
we come through to certain equations analogous to (8) and expressing the 
fundamental conjugate character of the solutions of (29) and (30), namely, 
the equations 


where D;,;,...:,;, 18 notationally the same determinant as that which is 
represented by this symbol in equation (8). 

Following still the method of § 3 it would be possible to name conditions 
on the coefficients a which are sufficient to insure that the first member of 
(31) shall have a value different from zero when p = oc. Since we have no 
need for the explicit form of these conditions we shall not take the space 
to derive them; we shall merely assume that we have before us a case in 
which the first member of (31) for p = a is always different from zero. 

It is obvious that we have particularly simple forms of the conjugacy 
conditions (31) for cases corresponding to those involved in equations (11), 
(12), (13); and that we may have modifications of all of the conjugacy 
conditions analogous to those indicated briefly for the finite case in the 
paragraph following equation (13). 

We are now led to consider the problem of the expansion of a function 
Zizi + i, Of r discrete arguments, each on the range 1, 2, 3, ---, in terms of the 
solutions of one set of systems, the expansion being of the form 


where the coefficients c; are independent of the variables 7;, ---, 7. If the 
function z has an expansion of this form and if certain changes of order of 


| 
| 
\ 
4 
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summation involved in the obvious process of reckoning are legitimate, we 
obtain readily for the coefficients c; the values 


r 

j h=1 


k= 2, 3, 


r 


= 


The discussion associated with equations (17) and (18) affords sugges- 
tions for the formation of special cases or for modifying the form of expansion 
in (32) and (33). 

Even in this transcendental case one may usefully employ a graphic 
representation of the sort described in § 5. One probably takes most natur- 
ally for an s-axis the real axis from 1 to infinity and erects perpendiculars 
at the points 1, 2, 3, --- on this axis, these perpendiculars to serve as the 
axes of coérdinates. One may then fill in the function by linear interpola- 
t tion as in § 5 and so make a function of a continuous variable. One might 
i conceivably carry out now a limiting process by which the number of 
perpendiculars in any finite stretch increases indefinitely and the greatest 
distance between two which are consecutive decreases to zero. One has a 
limiting case which may be employed to afford suggestions for dealing with 
differential equations of infinite order. Various modifications and exten- 
sions of this graphic representation are suggested by the foregoing treatment 
for the finite case, particularly that in § 7; but we shall not now pursue the 
matter further. 

Besides the generalization of the problem associated with (1) and (2) 
to that associated with (29) and (30) one might treat also that in which the 
number of parameters \ becomes infinite while the number of x’s and hence 
of y’s remains finite or that in which the number of )’s and the number of 
x’s and y’s simultaneously become infinite; but it appears unlikely that 
such generalizations will be as useful as that which is indicated in the fore- 


going paragraphs. 
A rather important problem, upon the resolution of which we shall not 
enter, may be suggested here, namely, the problem of the relative distribu- 
; tion of the zeros of linearly interpolated solutions of infinite systems. 
9. Variation of Parameters and Green’s Functions for Algebraic Systems. 
( —Let us consider a system of n non-homogeneous linear algebraic equations 
n+h 


(34) ait; = 1=1,2,---,n, 


and the corresponding homogeneous system 


n+h 


(35) == @, gm 1,2, 2; 
j=l 


| 
(33) ™ | 
| 
| | 
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each in n+ h unknown quantities, h being a positive integer, and let us 
suppose that the rank of the matrix || a;;|! is n. 

In what follows next we shall assume that the notation is chosen so 
that the determinant of the square matrix made up of the first n columns 
of ||a;;|| has a value different from zero. If the equations in each system 
are taken in a suitable order (and we shall suppose them to be so written 
already) it is obviously possible to combine them into a new equivalent 
system in the same unknown quantities in such way that the new coeffi- 
cients a;; have the value zero when 7 <7. We shall now suppose further 
that the original, and hence the new, matrix ||q@;;|| has the property that 
the determinants of order 1, 2, 3, ---, n in its lower right-hand corner are 
all different from zero in value. Then it is possible to make further com- 
binations of equations in either system so as to replace the system by an 
equivalent system in the same unknown quantities and in such wise that 
when all this is done the systems (34) and (35) take the forms 


h r 
(36) = Bi, = 0, 1,2, ++-, 2, 

r= 
where a;, ; and a;, ;,, are different from zero for every 7. For the equations 
in this normal form we have a ready method of obtaining the solutions of 
the non-homogeneous system in terms of a fundamental system of solutions 
of the homogeneous system analogous to, and indeed abstractly identical 
with, the classic method of variation of parameters in the theory of differen- 
tial equations, the latter being in fact a limiting case of the former. This 
method for the algebraic systems we shall now develop in a brief treatment. 

Let us denote by 


a fundamental system of h linearly independent solutions of the homo- 
geneous system in (36), and hence of that in (35). Then we seek a solution 
of the non-homogeneous system in (36), and hence of that in (34), in the 


form 


(37) 


(2 (h)(h 
{= 0,1,---,4— 1, 1,2, n+ 1, 


where the functions c? are to be determined. The equations in the system 


[eth — = 0, 0, 2, 1, 2, +++, 
(38) 
k=1 Qi, ith 
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are now obtained as follows: those in the first line come from (37) by 
comparison of that equation for a given value 7 + 1 of 7 with the equation 
gotten from (37) when ¢ is replaced by ¢+ 1; those in the last line are 
obtained by substituting the value of u from (37) into the non-homogeneous 
system of (36), the value of w+; for t = 0, 1, ---, kh — 1 being written as in 
(37) while that of w:;;, is taken in the form 


h 
k=] 


We look upon equations (38) as furnishing a system for the determination 
of the functions c{. 

The condition that (38), for every value 7 of the set 1, 2, ---, n and for 
every 8;, shall serve for the determination of the quantities in the square 


brackets is that the determinant A,, 


1 (2) (2) 
| 
| (2) (h) 
A; = | Vite 


(2) (h) 


shall be different from zero for every 7 of the set 1, 2, ---, n. But in § 6 
we saw that these inequalities are equivalent to the conditions that the 
determinant of every matrix of order n obtained from | |a;;|| by striking 
out h consecutive columns shall be different from zero. These are then 
necessary conditions for it to be true that every system (34) may be solved 
by the given method of variation of parameters in terms of a fundamental 
system of solutions of (35). These conditions, together with those already 
named in connection with the reduction to the normal systems (36), are 
also sufficient conditions for it to be true that every system (34) may be 
solved by the given method, as we shall now show by actually effecting 
the solution under the named hypoitieses. 

If we denote by AY the cofactor of the element in the kth column and 
last row of A; we have from (38) the relations 


+= 1,2, ---, 2. 


If we assign to c”, or to c“),, an arbitrary value c™, the equations (39) 
serve by recursion to determine uniquely the c? in terms of these arbitrary 


for fixed k and varying 7 is the arbitrary element 


values; and in fact c$ 
c plus a determinate function of i depending only on the b;, the a;; and 
the solutions of the homogeneous system. On substituting these values 
of the c’s into (37) we have the general solution u; of the non-homogeneous 


system in (36), and hence of system (34). 


| | 
1 
° ° ° | 
A(h) 
Qiy ith A; 
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This solution of (34) may obviously be written in the customary form 
(40) uz = CPAP + + + + 1=1,2,---,n+h, 


where 4; is a particular solution of (34) and the ce, ---, ce are arbitrary 
constants, this form of solution being of course always realized under the 
hypotheses of the first paragraph of the section. 

If we take systems of the form (36) where 7 runs over the infinite set 
1, 2,3, --+ instead of the finite set 1, 2, ---, n, it is easy to see that we have 
a process in every respect identical with the foregoing except that the re- 
currence relations which replace (38) are now to be solved by taking the 
first cf? only (since there is no last), namely c{”, as arbitrarily defined. 
Such a system is intimately connected with the theory of difference equa- 
tions. 

Let us now consider the question of adjoining linear conditions to the 
systems (34) and (35), these new conditions to be thought of as being 
analogous to the boundary conditions associated with differential equations 
for instance. Let these new conditions be taken in the forms 


n+h n+h 
(41) 
j=1 


j=l 
The question is on the possibility of choice of the ec in (40) and of the é” 
in the relation 
so that systems (41) shall be satisfied; that is, it is on the possibility of the 
systems 
h n+h n+h h n+h 
(42) Dice aya? + aii; = aja = 0, 
t=nt+1,---,n+s. 


In case s < h the latter system always has solutions ¢. It may also have 
solutions when s = h. When it has solutions the former system may or 
may not have a solution; its having a solution depends on the character of 
the a; and hence on the character of the b; in (34) and also on the character 
of the new a,; and 6; in (41). If s = h and if the last system in (42) has no 
solution except that for which the é“ are all zero then the first system in 
(42) has a unique solution. This last algebraic theorem (essentially obvious 
in its character) has for an analogue a fundamental theorem in the theory 
of differential equations (one which is not obvious in character).* 
Let us now consider the system 


j=l 


* Bocher, “Les Méthodes de Sturm,” pp. 19ff. 


| 

| 
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where the determinant A of the coefficients, A = |a;;|, is different from zero. 
Then the customary solution (which is unique) may be written in the form 


(44) =D Gic, 
i=l 


where 


A;; being the cofactor of the element in the 7th row and jth column of A. 
The function G;; has been called the Green’s function* for the homogeneous 
system 2a;;v; = 0 corresponding to system (43). For fixed 2 it satisfies 
every equation of this system except the ith, the first member of the 7th 
equation having for this value v; = A;;/A the value 1 instead of zero. Thus, 
by means of the homogeneous system, which is incompatible, we can define 
a Green’s function G;; of such sort that the solution of (43) for every set ¢; 
is expressible in the form (44). If the homogeneous system is not in- 
compatible it is easy to see that no function G;; exists having with respect 
to (43) the named property. 

If H;; is the Green’s function of the system adjoint to system (43) in the 
sense of § 2 it is obvious that H;; = G,;. 

This obvious matter for algebraic systems is mentioned here for its 
use in suggesting the requisite procedure for analogous properties of certain 
functional equations to be treated in later papers and especially for its 
use in transforming the expansions of § 4 into the form of contour integrals 
in § 10 immediately following. 

10. Condensation of Algebraic Expansions into Contour Integrals. Let 
us consider systems (10) with a single parameter \ and suppose that the 
determinant A(A) of the first system has n distinct simple roots AY, A”, 
--+,™ and that to each root corresponds a unque solution of the system 
(except for a constant factor). Denote by A;;(A) the cofactor of the element 
in the 2th row and jth column of A(A). Then the Green’s functions G;;(A) 
and H;;(\) for the two systems in (10) are respectively 

Ai;(A) Aji(\) 
Gij(A) = BQ)’ Hij(X) = Gyi(A) = 
These Green’s functions are analytic in the complex variable \ except at 
the zeros of A(A), that is, except for characteristic values \‘” of systems (10). 
Now we may write the Green’s function G;;(A) in the forms 
Gis) = 
* Bocher, Annals of Mathematics, (2) 13 (1911), pp. 71ff. 


) + Srij(A), k = 2, 


| 
n 
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where S;:;(A) is analytic in \ at X\=A™ and Ry; is the residue at 
4 = A of G;(A) considered as a function of the complex variable X. 
We propose to evaluate the residue R;,;; in terms of the solutions of systems 
(10). 

In the first place, we have - 
Reig = 
where A’(\“) denotes the value at \ = \™ of the derivative of A(A) with 
respect tod. The quantities 7; = R,,;, for fixed k and o, afford a solution 
of the first system in (10) forA = \. Likewise, the quantities y; = Rip, 
for fixed k and p, afford a solution of the second system in (10) for \ = \‘”* 
But all solutions of these systems for \ = \“ are constant multiples of 2? 
and y respectively. Hence we have 


(hk), (hk) 


where y; is independent of i and j._ It remains to evaluate yx. 
Now we have 
lim {A — — = 0; 


a=at*) 


whence it follows that 


(45) lim A—AM) DD — = 0. 
} i-1 j=l 


A=A‘) i=l j=l 
But 


n 


(cas + = — AM) Dua, 


so that from the fundamental property of the Green’s function we have 


Putting this value of x? into (45) we have readily 


1 — ye = 0 
4#=1 j=l 


With the value of ;, afforded by this relation we have for the residue Rx;; 
of Gi;(A) at \ = A in the complex plane the value 


a (7) 


Vi 
y® 
p=1 o=l 


From this result we have for the kth term in the expansion of z; in terms 


n 
n n 
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of x! by (16) and (17) the value 


i=l j= k 

jie 

j=l 
where I’; is a contour about \“ inclosing no other characteristic value of 
system (10). Then, if I is a contour inclosing all the characteristic values 
of systems (10), we have 


pi=i j=1 
Thus we replace the expansion for z; by a contour integral. It is clear that 
any set of terms of the expansion may likewise be replaced by the same 
integral taken about a contour inclosing appropriate corresponding zeros 
of A(A). 

Contour integrals of this character were first employed by Birkhoff.* 
They served admirably in the convergence proofs associated with the 
expansions of functions in terms of functions defined by differential systems 
these systems being of such character as to be realized as limiting cases of 
the algebraic systems here in consideration. We shall find of fundamental 
value in our general investigation various limiting forms of the contour 
integrals developed in this section. 

It is possible (but not important for our purposes) to modify and extend 
the foregoing analysis so as to avoid the restriction that the zeros of A(\) 
shall be distinct. 

Let us now pass to the more general systems (1) and (2). For any 
given fixed value of h let A; (Ai, A2, ---, Ar) denote the determinant in the 
first member of (3) and let Ajj; (Ai, Ae, «++, Ar) denote the cofactor of the 
element in the 7th row and the jth column of this determinant. Then the 
Green’s function G,;; of the hth system in (1) has the form 


Anij(M1, de, dr) 


Cris An(A1, » 92 as 


We assume for the present purpose that systems (1) [and also systems (2) | 
are so restricted that the solution of each for characteristic values of the \’s 
is unique (except for a constant factor) and that 


* Transactions of the American Mathematical Society, 9 (1908): 373-395. 
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exists and is finite, where dA”, ---, A” is any set of characteristic 


systems (1) and (2). Then the function 


i k) k) k 


of the complex variable ,, has a pole of the first order at the point \, = X”. 


Let its residue at this point be denoted by Ryai;. 


In order to evaluate the residue Rxx:; we apply the earlier result of this 
section for R;:; to the system which results from the hth system in (1) on 
replacing each A; except A, by the corresponding characteristic value \%. 


Thus we have 


k), (k 
= 


where x, is independent of i and 7. Then we have 
k) , (k) 
Tl Rasy, = vel] 
h=1 h=1 


where y;, is independent of the 7 and the 7. 
Then if we write 


(A1, de, Ar) = II Gni,j, 15 da, 


we have 


k k 


r 


whence it follows that 


Ar) 


t=) 


Employing the Ath system in (1) and proceeding by the method used 


h=1 


values for 


Ak) (k) . 


ay}, 


A=1 


II 


hins 


earlier in this section for evaluating x in terms of the Green’s function, 


we have 
s=1 t,=1 


From this we may form the product 


Ak) k 


and simplify its expression by taking the limit for any convenient approach 


9 i 
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of Ar to the characteristic values ---, A“. Then, in a certain 
range of cases (if not always) it turns out that the value of the product in 
question is that of the first member of equation (48). We confine attention 
to the cases in which this new condition is realized. Then we have 


ny ny 


Thence it follows that 


r 
(ke), (ke) 
II 


From this result we have for the kth term of the expansion of 24... 
afforded by (14) the value 


ni n,. nN, 


where T',; (h = 1, 2, ---, r) is a contour in the d;-plane about the point 
\ and containing in its interior no other characteristic value of \,. If 
we replace the contour ',; (h = 1, 2, ---, 7) by Ty, a contour which includes 
within it all the characteristic values of \,;, and perform the same multiple 
integration about such contours, we shall have the value of the function 
Znty..t,. It is clear that we may form similarly the contour integral for 
any given partial sum of the series for 2;,...., afforded by relation (14). 


UnIvERsITY OF ILLINOIs, 
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ALGEBRAIC THEORY OF THE EXPRESSIBILITY OF CUBIC FORMS 
AS DETERMINANTS, WITH APPLICATION TO 
DIOPHANTINE ANALYSIS. 


By L. E. Dickson. 


1. It was proved geometrically by H. Schréter* and more simply by 
L. Cremonaf that a sufficiently general cubic surface f = 0 is the locus of 
the intersections of corresponding planes of three projective bundles of 
planes: 


+ Alyy + = 0, + Aloo + = 0, kls1 + + = 0, 


where x, A, w are parameters and the /;; are linear homogeneous functions 


of 21, --+, x4. Hence the surface is expressible in determinantal form 
\l:;| = 0. Solving the three equations, we get 
= fi fo ifs 


where f; is a homogeneous cubic function of x, A, nu. We thus secure a para- 
metric representation of the points of the surface, which is therefore uni- 
cursal. 

In case all of the coefficients of the /;; are rational, we have the complete 
solution in rational numbers of the Diophantine equation f = 0. However, 
there exist cubic equations f = 0 whose rational solutions involve three 
parameters homogeneously such that f is not expressible rationally in 
determinantal form (§ 12). 

With the application to Diophantine analysis in mind, I here discuss 
algebraically the problem to express a given cubic form f as a determinant 
\1;;! of the third order. Since we are interested ultimately in the case in 
which the coefficients of the /;; are rational and since three linear functions 
of four variables vanish for values not all zero of the variables, we shall 
assume that f vanishes at a known rational point. Then the coefficients 
of the /;; are shown to be expressible rationally in terms of a root of an alge- 
braic equation whose leading coefficient is not zero if f = 0 has no singular 
point. The existence of a rational root may be decided by a finite number 


of trials. 
The corresponding binary problem is solved by the identity 


* Journal fiir Mathematik, Vol. 62 (1863), p. 265. 
t Ibid., Vol. 68 (1868), p. 79. Cf. Clebsch, ibid., 65 (1866), p. 359. 
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ax+by y 0) 
ax? + bx’y + cry? + dy = —cy x ~ 
dy 0 2 | 


In the ternary case we obtain simple geometrical criteria that a cubic 
curve with a rational point be expressible rationally in determinantal form. 
There exist determinants of the third order which vanish for no rational 
point, for example, the determinant A(£) of the general number 
€= a+ y0+ 2¢ of the algebra having the units 1, 0, ¢ = 6, where 0 is a 
root of a cubic with no rational root. But ternary forms without rational 
solutions have little interest in Diophantine analysis and will not be treated 
here. 

Given one representation of a form as a determinant, we may derive an 
infinitude of representations by the familiar operations which leave a 
determinant unaltered in value. The more definitive problem treated here 
is to find a representative of each class of equivalent matrices with linear 
elements which have a given determinant. This raises the question of the 
number of classes of such matrices. 

In another paper* I prove that every binary form, every ternary form, 
and every quaternary quadratic form are expressible in determinantal form; 
while, apart from these and the quaternary cubic, no further general form 
has this property. The rarety of such forms justifies the present investiga- 
tion of cubic forms with attention to rationality. The theory covers every 
cubic form, not merely a sufficiently general one. 

2. The given rational point on the locus can be transformed rationally 
into (1, 0, ---, 0). Then the locus is f = 0, where f = 2°f1 + af2 + fs, 
where f; is of degree j in y, z, ---. If f1 is not identically zero, we take it 
as the new variable y. Then by adding to 2 a suitable linear function of 
y, 2, +++, we may delete rationally from f: all terms with the factor y. 
Hence 


(1) f= xy t+ afet+ fs, 


where the quadratic f2 is free of 2 and y, and the cubic f; is free of x. But 
if f; = 0 and f. = 0, we may transform f2 into ay” + q, where a + 0 and q 
lacks y. By adding to 2 a suitable linear function of y, 2, ---, we may 
delete rationally from fs; the terms with the factor y? and obtain 
f = + yr+s, where q, r, s lack both and y. Taking az asa 
new x, we havea = 1. Interchanging 2 and y, we obtain a form of type (1). 

3. Let the form (1) equal to a determinant of the third order whose 
elements are linear forms with rational coefficients. We may assume that x 
occurs in the first element of the first row and that its coefficient is unity. 


-* Transactions Amer. Math. Soc., April, 1921. 
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The remaining elements in the first row and first column may be assumed 
to be free of x. By interchanges of the last two rows and last two columns, 
we may assume that the second element of the second row contains 2; 
its coefficient may be made equal to unity. Hence we may take 

+ ls ls 
(2) f=| 

| l; ls 
where each /; lacks x. After subtracting multiples of the last row and last 
column from the remaining rows and columns, respectively, we may assume 
that /3, Js, lz, lg lack 2 and y. Thus by the terms linear in 2, 


(3) + 0, = = Isls = fo. 


The interchange of the first two columns and the first two rows of (2) 
corresponds to the substitution 


S = (lola) (lal). 
The interchange of rows with columns corresponds to 
= (Isls) (Isl7) (Ils). 


Add the products of the elements of the second row by k to the elements 
of the first row and then subtract the products of the elements of the 
first column by k from the elements of the second column; we obtain a 
determinant of the same form with the same J,, J, 17 and with 


li +. = ls kl, -+ kl; — k?l,, l; = l; 


(4) 


Proceeding as before with the words first and second interchanged, we 
obtain a like determinant with the same 2, /3, 1s and with 
= 1, — kh, 14+ kl, — — Ila, = 
The last result may be obtained by transforming (4) by S. 


4. For 3 variables (1) becomes 


(6) C= apt Byst be. 


The coefficient of x* in the Hessian of (6) is — 8e. Hence (1, 0, 0) is a 
point of inflexion of f = 0 if and only if e = 0. If e + 0, we multiply y 
by e? and 2 by 1/e and obtain a form of type (6) with e = 1. We seek the 
conditions under which (6) can be expressed as a determinant (2), where 
ls, Ig, lz, lg are multiples of z. 


| 
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First, let either 1; = 15 = 0 or 1; = 1; = 0. Then (2) is the product of y 
by a quadratic which evidently vanishes when x = — i, 1, = 0. Hence 
(6) must be the product of y by a ternary quadratic g which vanishes at a 
rational point. Conversely, the product of such a q by a rational linear 


form / can be transformed* rationally into one of 


Kay —k22))=|0 l(ay? + byz + =|0 ayt+bz 2z 
O kz y (0 — cs yl 


However, the product of a linear form, say y, by an arbitrary ternary 
quadratic form Q is expressible rationally as a determinant. If Q lacks*a?,* 
it vanishes at (1, 0, 0). In the contrary case we may delete the terms in 
ay and xz by a transformation on x and use identity (7). 

Second, let /3, Js, 17, ls be not all zero identically. Since S, 7 and ST 


replace /; by l/s, 13, ls, we may take /; = 0. Hence we may set 1; = z. By 


(4), we may take 43 = 0. By (3), — 132 = e2?, whence J; = — ez. 
We first treat the case e = 0. Then /; = ls = 0, 1; = z. Hence shall 
dy l, O | 
C = | ls ls ly 
| y | 


where /; is a multiple bz of z. If l, = ay, C = y(— [| — ayly + abz?), which 
requires 6=0 m C and is then satisfied if a=1, —ay—z, 


so that we have the identityt 


y 0 | 
(7) y(a? + ay? + Byz + =|—ay— Bz x 
| z O y | 


Next, let /2 contain z, whose coefficient may be divided out of the second 
column of determinant C and multiplied into its second row. Then 
l,=2z-+cy. By (5), we may delete z from J, and set 1; = ay. Write 
ly = gy + hz, i, = tz. Then shall 
ay z+tcy 0| 
C=|gythz —ay tz| 
z 0) y | 
= — (a? + cg)y® — (g + ch)y’z + (ct — + te’, 
whence 
= 6, h = — y, g = —ch— a+ecg=—a. 


The last relation furnishes the condition 


+ Dickson, “ Algebraic Invariants,” 1914, p. 24. 
7 If (6) with e = 0 has any linear factor, it has the factor y. 
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(8) ye + = 


which is solvable in rational numbers if and only if (6) with e = 0 contains 
the rational point* (a, 1, — c). The latter may be any point not on the 
inflexion tangent y = 0; if 6 + 0, it is any point except (1, 0, 0). 


Next, let e = 1, whence 13 = — z. We have i = kz. By (5), we may 
take = 0. Thus (2) becomes 
(9) a—l 0 | = gy + a2? — y(li? + bela) — 
| 


z 0 y 


Identify this with (6). By the terms in 2°, we see that /; = ay — 8z. 
Set lo = cy + dz, ly = gy + hz. Then the conditions are 


cg= —-a-— a’, dg + ch = 2a6 — 8B, dh = —y— & — a. 
Also, dg — ch = 2r must be rational. From 
(dg — ch)? = (dg + ch)? — 4cg-dh, 
we obtain the condition that 
(10) — 4a* — 4a°y — 4a(a + BS) + B 4a(y + = 


shall have rational solutions a, r. Then we have rational values of dg = A, 
ch = B, dh = C, cg = D such that AB= CD. If A +0 or C +0, we 
may take d= 1, g= A, h=C, c= D/A or B/C. If B+0 or D+0, 
we may take c = 1,g = D,h = B,d = C/B or A/D. 

To interpret (10), replace 2 by /; = ay — 6z in (6) with e = 1, and re- 
move the factor y. We obtain a quadratic in y : z with a rational rootf 
if and only if (10) has rational solutions. Then (6) has three rational 
points of intersection with x = ay — 6z for some rational value of a, i.e., 
for some line other than y = 0 through (— 6, 0,1). The latter point is the 
tangential of P = (1, 0, 0), i.e., the new point in which the tangent y = 0 
at P to the cubic curve meets the curve. 

THEOREM. Any reducible cubic curve is expressible rationally in deter- 
minantal form. <An irreducible cubic curve with a rational inflexion point is 
expressible rationally in determinantal form if and only if it contains a further 
rational point. A cubic curve with a rational point P not an inflexion is 
expressible rationally in determinantal form if and only if it has three rational 
points of intersection with some line, other than the tangent at P, through the 
tangential of P. 

5. For four variables, we may set fo = az?-+ bw? in (1), where b = 0 


* For which z = l,, lp = 0, whence the elements of the second column of (2) are all zero. 
t Which makes /: = 0 in accord with the second column of (9). 


| 
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ifa=0. First, let /3, ls, 17, ls be not all free of z, the contrary case being 
treated in $9. In view of substitutions S and T of § 3, we may assume that 
I, contains z; its coefficient may be made unity by removal of factors from 
the last row and last column of (2). By (4), we may assume that /s = ww. 
We here take » + 0 (treating » = 0 in $10). By dividing the elements 
of the second column by v and multiplying those of the second row by 2, 
we have /s = w. Then by (5) we may take /; = z. Hence,* by (3), 


ls = — az+ gu, le = — gz — bu, 
ath ls ls ls | 
w y wy 
Set 


(12) h=cytd+t+ew, k= my+nzt+ pu, 
(13) fs = Ay® + By’z + + + Eyzw + Fyw? + 


+ Hz2?w + Jew? + 
Then (11) is identical with (1) if and only if 


—c’—hm= A, — 2cd—hn—jm=B, — 2ce—hp—km=C, 
— @—jn—gh—ac=D, — 2de— jp — kn+ 2gc — am— bh= E, 
(14) —e—kp+be+ gm= F, —ad—g=G, 
2gd — an — bj — gk — ae = H, 
bd + gn + 2ge — ap — bk = J, be + gp = L. 


We shall designate these equations as (A), ---, (L). 

First, lta = b=0. Ifg = 0, then /; = J, = 0 and (11) is the product 
of y by its minor. Postponing this case to the end of § 11, we have g + 0. 
Let dg = 6, eg = «. Then equations (L), ---, (D) give 


gp=L, gj=—G, gn=JS—2¢, gk=26—-H, 
PE—LG—HJ+ 2eH+ 26) — 26e. 


Inserting these values into the product of (B) and (C) by g’*, we see that 
the terms of the third degree in 6 and ¢ cancel, giving 


(15) + 26H + 326 + — 3LG — HJ) + €(2g2D — 
+ gB — g(FG + DJ) + GJ?+ GHL = 0, 
(16) 38L + + @H + — 4HL) + — 3GL — HJ) 
+ gC — g(FH + DL) + HL + GJL = 0. 


*It is proved in $13 that no further normalization of our determinant is possible 
and that the two excluded cases are truly exceptional. 


| 
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Multiplying (A) by 39°, we obtain an equation whose terms of highest 
degree are — 66°L — 66°eJ — 66H — 6G. Hence by adding to it the 
product of (15) by 2€ and the product of (16) by 26, we obtain the quadratic 


F —5HL+3J*) + D—5GI+ 3H?) + 
+ 6[29'C+ 9°(3EJ 
+ g°(3EH—2DJ—8FG)+5GHL+ 


(17) 


The true resultant of three ternary quadratic forms wu, v, w in 2, y, 2 is 
known* to be the determinant of the coefficients of 2”, ry, xz, y’, yz, 27 in 
U,V, Jz Where 7 is the Jacobian (functional determinant) of w, 2, w, 
while j, is its partial derivative with respect to x. The resultant is of the 
fourth degree in the coefficients of each form u,v, w. To find the degree in g 
of the resultant R of our forms (15)—(17), we shall determine the terms of R 
of highest degree in g. Writing the functions (15)-(17) homogeneously in 
5, e, 7 and retaining in each coefficient only the term of highest degree in g, 
we have 


+ 2eJ + @PH+ erPE+ 7’9'C, 
G(OF+ bE + 2irg?C + 2erg?B+ 


Thus R = (g*)*r, where r is the resultant when the factor g? is omitted 
from the last form. We make the transformation of variables 6 = w, 
€ = 2, gr = y and obtain the partial derivatives of f3, given by (13), with 
respect to 2, w, y, respectively. Hence the resultant of the final forms is 
the discriminant A of f;. Let J denote the Jacobian of these derivatives, 
and 7 denote the Jacobian of our forms having the resultant r._ By a general 
theorem, j equals the product of J by the Jacobian g’ of w, z, y with respect 
to 6, Now J is independent of Hence the exponent of = in 
any coefficient of 7 exceeds by unity the exponent of 7 in the term. ‘This 
is therefore true also of the derivatives j;, j., while in 7, the exponent of M 
in any term exceeds by 2 the exponent of 7. We remove the factors M, 
M, M? from the last three functions. We now have six quadratic forms in 


5, e, T, such that M does not occur as a factor of a coefficient of 5°, de, &’, 
while M (but not M?) is a factor of the coefficients of 67 and ev, and MW? 
(but not J*) is a factor of the coefficient of 77. Their determinant thus 
has the factor M*, but not M°®. Hence r= M*-M*A. But R= M’*r. 
Hence R= MA. Thus g? is a root of an equation of degree 12 whose 
leading coefficient is the discriminant of f3. 

” Salmon, Algebra, § 90. 
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The constant term of this equation is not identically zero. It is in fact 
not zero when J = H=0, LG +0. For, then equations (15)-(17) 
become for g = 0 

3G(e — 6L) = 0, 3L(8& — eG) = 0, 3LG(— + = 0 


By the first two, dé = 0 or LG. Hence the resultant is not zero. 

The coefficients in (15)—(17) are independent of g if and only if A, B, 
C, D, FE, F are all zero. If they are zero the resultant is free of g and, as 
just shown, is not zero if J = H = 0, L@ + 0. 

THEOREM. The problem to express xy + f3 rationally in determinantal 
form, where fs is a cubic form in y, 2, w, depends completely upon the deter- 
mination of a rational square which satisfies an equation of degree 12 whose 
leading coefficient is the discriminant of f3 and whose constant term is not zero 
identically. The equation reduces to its constant term uf and only if fs lacks y. 
In particular, if LG + 0, ay + Gz? + Lw? is not expressible in determinantal 
form under our present assumptions (cf. § 10). 

The equation for g involves only even powers since this is true of (15)- 
(17). A more fundamental reason is given in § 14, where there occurs an 
example in which the equation for g’ is a cubic, the discriminant of fs 


being zero. 


To give another example, take /, = w, = — y, 4 =0,g= 1. Then 
x u 
w — = w+ wy + 
w y 


with no rational singular point. Since C = D = L = 1 and the remaining 
coefficients of f3 are zero, (15)—(17) become 
= 0, 30° = — + (2g* — = 0. 


We desire that g + 0. The only real solution is 6 = e = 0,g? = 1. Hence 
the only rational representations as a determinant of type (11) are the 
above determinant and that obtained by changing the signs of the elements 
other than 2 in the first two rows. 

6. We now discard the assumption that a = b = 0, but assume that 
g+0, A= g+ab +0. It will prove convenient to introduce d; = Ad, 
e; = Ae in place of dande. The last four equations (14) give 


gAp = AL — bey, gAj = — AG — ad,, 
gAk = 2gd, + ae, + 6, gAn = bd, — 2ge, + «, 


where 


§= — PH + g(bG — aJ) — e= gJ + g(al — bH) + BG. 


| 
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In (D), (E), (F), the determinant of the coefficients of h, m,cis 2gA. Hence 


2g4h = aN-+ Bag — ya’, a=-D-d—jn, N = ab + 24’, 
2gAm = — ab? + Bbg + yN, = — E — 2de — jp — kn, 
2gAc = Fret+kp. 


Replacing p, j, k, n by their values, we get 
= di(ab — — 2dyeyag + di(AbG + ae) — 2e,AGg + AGe — 
= — 2dibg + 2dye:(g? — ab) + 2ejag + d,(ALa — bb — 2ge) 

+ — ae — AGb) — de — A*Eg? + A°LG, 
g’A?y = ei(g? — ab) — 2dyeybg + 2d,ALg + e:(ALa — bd) + ALS + 


Mulitplying equation (B) by 29*A* and inserting the preceding values of 
h, m, c, n, j, we obtain 


(18) 2djAab? — 2dieAabg + — 
+ dj{A(5Gb?9? + 3Gab? — 2Lag? — 4La*bg + 3abe) + 29°(b5 + ge)} 
+ 2de1{(ALa? — 6N)(g? — ab) — 2AGbgN — 2a*bge} 
+ ej {2ALa'g + (AGN + ea?)(ab + 3g?) — 26gaN} 
+ di{— 3A°LGgN — ALa?(ge + 2b5) + 2A?Eg>’ — 2A?Fa*bg? 
— 2A°Dbg?N + + 4AGbeN + + 
+ (29° — abg)de} 
+ e1{ (ade — A?LGa)(ab + 49) + LAa?(295 — ae) + 
+ + 2A*DgiN — A°G?b?g — 4AGgeN + AGab?s — a°ge*} 
+ 2A‘g*B + (A?Gb? + AeN) (Ge — ADg’) 
+ (age — AGbg)(A°LG — A?Eg? — 
— (A°GN + Aa%e)(Lé + AF9?) = 0. 


Similarly, from (C) we obtain an equation 


(19) — 2diAb’g — 2die,Aab? — 2dye{Aabg — 2e{Aa*b + ---, 
which may be deduced from (18) by the substitution* 
(20) (BC) (DF) (GL) (HJ) (ab) (di, — €1)(9, — 9), 


which leaves A, LE, A, N, 6 unaltered and induces 


* (d, — e) replaces d by — e ande by — d. 


i 
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(6, — e)(a, — y)(c, — ¢)(hm)(kn) (jp). 


These substitutions and (zw) interchange J, with J,, 1; with Js, and change 
the sign of /; and hence interchange the first two rows and first two columns 


of the initial determinant (11). 
There remains the first equation (A) in (14). Its product by 49°A5 


may be written in the form 

(21) 4g°A°A — g'A*(ab — ya)? + + = 0. 

Its terms of the highest degree in d, e; are (3g? — ab)A, where 
= A(dib? + 2dieiab + eja?). 


The terms of the fourth degree in bd,(18) — ae\(19) are 2abdA; those in 
gei(18) + gd,(19) are — 2g’\. Hence by adding multiples of (18) and (19) 
to (21), we may cancel the terms of the fourth degree. Hence we have 
three cubic functions of d,, e;, whose true resultant R is known to be express- 


ible as a determinant. 
We proceed to find the terms of R of maximum degree in g. From the 


aggregate of terms in (18) which multiply each djei we omit all terms not of 
the highest degree in g, remove the common factor 2 and obtain 


(22) dig’ab? — dieg’ab + — ejg?a? + dig’J + 
+ + dig®E + 2e19°D + 
Applying to this the substitution (20), we obtain 
(23) — dig%b? — die.g’ab? — — + 3dig®L + 
+ + + + 
The resultant of these two functions and a third function 7 of d, e; is 
the same as the resultant of their sum and difference and 7. Hence we 


may omit the terms involving only g?. Remove the factors g*, make the 
functions homogeneous in d,, €;, 7 and set g*r = ¢; we get the forms free of g: 


(24) — die,ab — + + 2dieidH + + di¢?E + 2¢1¢°D + 
(25) — dib? — diab + 3dibL + + + + + PC. 
In (21) the terms of maximum degrees in g are 
(26) 3dig'b? + 6dieigtab + 3etgta? — SdigiL — Sdieig'J — Sdyeig'H 
— S8eig'G — 4dig’F — — 4eig’D + dig'(4JE — SLD) 
+ (4HE — 8GF) + 


The terms in g'! near the end do not contribute to the part of the resultant 


i 
| 
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of highest degree in g, since by omitting them (or replacing them by g'® 
with the coefficient zero), we shall obtain a resultant not identically zero. 
Also divide by 4g’, make homogeneous in dj, ¢;, 7, and set g’r = ¢ as before; 
we get the form free of g: \ 
(27) 2dib? + 3dieiab + — 2digL — — — 
— di¢’?F — — ej¢’D + 
Now (24) and (25) are the partial derivatives with respect to e; and d, of 
(28) A¢' + Cd,¢° Dei¢? -+- Ed,e,¢? Fdi¢? +. Geid Hdyeid 
+ Jdied + Ldi¢ — 4dib? — Sdiejab — 
= $f3(¢, e1, ci) — 4(dib + eja)’, 
where f3(y, 2, w) is given by (13). To (27) add the product of (24) by e; and 
the product of (25) by d;; we get (28). Hence the resultant of (24), (25), 
(27) equals the resultant of (28) and its partial derivatives, i.e., the dis- 
criminant of (28). This discriminant is not zero for all values of A, ---, 
L, a, b. In fact, the general quartic curve becomes 2K — 2°y? = 0 when 
referred to a triangle of reference whose side z = 0 is a bitangent and whose 
sides 2 = 0 and y = 0 are any lines through its two points of contact. 
Evidently (28) is of this form if x, y are the factors of 3(djb + eja). 
The quartic curve (28) is a plane section* of the tangent cone to the 
cubic surface (1), viz., 


(1’) + x(aei + bd) + e1, di) = 0 


whose vertex is ¢ = d; = e, = 0,2 = 1. A point on (1’) will be a singular 
point if the partial derivatives with respect to 2, e1, d, vanish: 


9 2 E 
+ aei+ bd?, 2axer + 2abd + 
1 1 


Let ¢ + 0 and substitute the value of x for which the first function vanishes 
into the others and (1’), and multiply each result by ¢; we get 
2 2 Ofs 2 2 Of 
(29) — ae;(ae, + + 950, bd; (ae, + + Pad,’ 
(ae; + bd)? + 


which equal functions (24), (25), (28), respectively. Since a general cubic 
surface is reducible to (1’) by § 2, and has no singular point, we again | 
conclude that the resultant of (24), (25), (27) is not identically zero. 


Miller, Blichfeldt, Dickson, Finite Groups, 1916, p. 352. 


. 
| 


Dickson: Expressibility of Cubic Forms as Determinants. 113 


Multiply the first function (29) by e; and the second by d;, add and apply 
Euler’s theorem Le,0f3/de, = 3f3; we get 


€1(24) + d,(25) = — + bdi)?+ (3 


Subtract this from 4 times e;(24) + di(25) + (27) = (28), given above. 
Hence 


3e1(24) + 3d1(25) + 4(27) = 


Hence we may replace (27) by the cubic function 0(@f3)/d¢. This may 
therefore be obtained by deleting the factor g* from 


3e19(22) + 3dig(23) + (26), 


after the three functions are abridged as above by omission of terms in g? 
and g''. We also deleted the factor g* from each of the abridged cubics 

22), (23). But the resultant is of the ninth degree in the coefficients of 
each form, and hence equals the product of (g*)°(g*)°(g*)® = g® by the re- 
sultant r of the three cubic forms homogeneous in d;, ¢;,7. Their Jacobian 
j equals the product of the Jabocian J of the equivalent forms in dj, e,, 
@ = g’t, by the Jacobian g’ of di, e:, @ with respect to d;, ¢1,7. But J is 
independent of g. Hence the exponent of M = g? in any coefficient in 7 
exceeds by unity the exponent of 7 in the term. The same is therefore 
true of the derivatives ja, and j.,, while in j, the exponent of W/ exceeds by 
2 the exponent of 7 in any term. We remove the factors @, MW, M? and 
obtain three quintic functions in which the exponent of M in any coefficient 
equals the exponent of t in the term. The same is true of each of our three 
cubic forms. If we multiply each of them by dj, dye, ej, dit, e:7, 7, we 
obtain 18 quintic forms. The determinant (of order 21) of the coefficients 
in these and the former three quintics is known* to be the resultant p of the 
three cubics. In each of the three equations obtained by use of the multi- 
plier 7°, the exponent of M in any term is 2 less than the exponent of 7; 
we multiply these equations by M?. Similarly we multiply by M each of 
the six equations obtained by use of the multipliers d)7, eyr. Now in all 
21 equations the exponent of / in any term equals the exponent of r. 
Hence in the new determinant the elements of any column contain the same 
power of M as a factor, viz., the exponent of 7 in the corresponding term 
der", 1+ m+n= 5. Hence the total exponent of the power of M 
which divides the determinant is 5 + 4(2) + 3(8) + 2(4) + 1(5) = 35. 
But our new determinant equals p(M?)?M®. Hence p is the product of 
M*/M” = M*® by a constant \ + 0. Thus, accounting for M* removed 


* Salmon, Algebra, § 90. 
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above, r = M?. Hence the resultant of our initial equations is of degree 
90 + 3-27 = 171 ing. 

THEOREM. When the quadratic form f, in 2, w does not have rational 
factors, the problem to express x’y + af2 + f3 rationally in determinantal form 
depends completely upon the solution of an equation of degree 171 whose 
leading coefficient 1s not zero if the surface has no singular point. 

Perhaps this equation is reducible since its degree is high in comparison 
with the degree of the equation upon which depends the determination of 
the 27 ruled lines on the surface. 

Material simplifications arise if G= H = J = L = 0, whence y is a 
factor of f3, and y = 0 cuts the surface in ruled lines. Since 6 = ¢ = 0, 
(18) and (19) have no quadratic terms. Dividing 6d,(18) — ae;(19) and 
€1(18) + d,(19) by 2, we get 


(30) Aab(dib + eja)? + — Fab? — Db?N) 
— ej A’g?(Eg’a + Fa?N + Da’b?) + — Fag + Dbg) 
+ A‘g*(d,Bb — e,Ca) = 0, 


(31) — Ag(d3b + efa)? + d2A%g3(— Egb + FN + Db?) 
+ ejA’g?(Ega + Fa? + DN) + 2d,e,A*g*(Eg + Fa — Db) 


+ Atg(d,C + eB) = 0. 
Dividing g(30) + ab(81) by A?9*, we get 


(32) (bdi — ae?)[E(g? — ab) + 2Fga — 2Dbq] 
+ 2dye,[2Eabg + (g? — ab)(Db — Fa)]+ A%bdi(gB + Ca) 
+ A*ae,(Bb — Cg) = 0. 
We obtain a second quadratic from [ (21) + (30) + 3g(31) ]/A’9?: 
(33) di 2Ebg*+ Fp+ D(ab*— g’b*) |+ —2Eag?+ F (a*b— g?a”)+ Dp | 
+ E(g?+3ab) — Bb) 
+ A’g'(E?—4F D+ 44 A) —A’g?(Db+ Fa)?=0, 
where p = a*b? + 5abg? + 29. Here the coefficients of dj and e{ will be 
proportional to b and — a, as in (32), if and only if aF + bD = 0, and then 


the same fact holds in (30) and (31). If also B = C = 0, no first degree 
terms in d;, e; occur in our equations. Then (32) and (33) give 


— ae{)Q=4Al[ Eabg— dieiQ=AL— E(g’?—ab)—4F ag ]SQ’, | 
Q= (9? — 3ab)t, t= S = 24A+ Ft. 


| 
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From the squares of these, we get (if* ¢ + 0) 
(bdi + ae{)?(g? — 3ab)Q = 
We obtain a second expression for the left member by use of 3(30) + (31): 


— + ae})*(g? — 3ab)A + — ac?) Ez +5 Gab + 20) | 


+ + 3abg) — 8Fag?] = 0. 


Multiplying this by Q, inserting the preceding values, and cancelling the 
common factor — 2A‘g*, we get either S = 0 or 


g'(4abA + t) + 2abg?(4abA — t) + a?b?(4abA + t) = 0, 


whence 


{g?(4abA + t) + ab(4abA — t)}? = — 16a°b*At. 


Hence the number of rational values of g? is 3 or 1 according as — abAt 
is or is not a rational square. If also d= D=—1, E=2, F=1, 
a=b6=1,then g=+1,+1 +42, 41 — v2, and the surface has no 
singular point. 

7. In §$ 7-11 we shall treat the special cases which were excluded in 
§§ 5-6. In each case the determinantal surface has a known ruled line 
whose equations are so simple that its occurrence on the given surface can 
be detected by inspection. When the given surface has the line as a ruling, 
its representation as a determinant is a much simpler problem than that 
treated in §§$ 5-6. Accordingly one should first ascertain whether or not 
the given surface falls under one of these special cases. 

Consider the case g?-+ ab = 0 which was excluded in §6. If a= 0, 
then b = 0 ($5) and g = 0, J; = Js = 0, so that (11) is the product of y 
by its minor, a case treated at the end of §11. Hence leta +0. We 
first treat the case b + 0. Then 


bw? =bZW, Z=wt W=w- lZ=gW, Ilk=—bW. 


Introduce Z and W as new variables in place of z and w’and let the new fs 
be given the same notation (13). We have therefore to consider a form (1) 
in which f. = bzw. Interchanging z and w if necessary, we may assume 
that, in (2), /; contains z. Proceeding exactly as in § 5, we may set /s = w, 
By (3), = tz, 13 = (1 But J; and must be dependent 


and hence one of them zero. Applying substitution S, we may take /, = 0. 


*If t = 0, either Sg? = 0 or E = F = 0, since the determinant of the coefficients of 
E and F equals — A?. For E = F = B =C =0, D = 0, (83) gives A = 0, and f has 
the factor z. 
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Thus we have (11) with 4, = 0, /13 = w. It remains to identify the final 
determinant in (11) with f;. Set 4 =cy+1, my+s. By the 
terms free of y, 


w(wr4 + = G2 + + Jew? + 
whence G = 0, A, = Lw+ nz, \; = Hz+ (J — n)w. From the remaining 
terms we remove the factor y and have 
— T= bly + w(wm + 2c) = Ay? + Byz + Cyw + D2 + Ezw + Fw. 
The resulting conditions (A), ---, (/) are 
—?—mh=A, —2cH—mj—nh=B, 
—-H’?—n=D, m—Lk— (J =F. 


First, let n + 0. From (D), (£), (F), (B), 


m= —+—, | 
J n n n>’ n n 


where 


N= D+ P= LN — n(E+ 2HJ) + 2Hn’, 
R= nvF+ (J — n)?n?+ LP. 
To clear the denominators of n multiply (A) by n®, (C) by n‘, and use the 
abbreviations 
N= D- BR’, S= E+ T= F+ J*+ 20H, 


R= 2n?J 4+ WT — nlS+ PN. 
Hence 


(34) — 2n*LH + n?L?N) + cL2R(nLN — — n*LB] 
+ — Rn’B+ n'A = 0, 
(385) en? L+ e(— n° + — LS + + — nS + 2LN) 
— = 0. 
Retaining in each coefficient only the highest power of n, we get 
— 2en'H + n8N = 0, en? L — + = 0, 
whose resultant is n*N. Hence n satisfies an equation of degree 24. 
Second, let n = 0, L + 0. This case occurs only if G = 0, D = — H’. 
Use the abbreviations S = E+ 2H/J, B= F+ J. Then (£), (F), (C) 


give 
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m= B+ Lk, Dh = — C — 2cJ — k(B + Lk). 


By (B), c=A/u, X= — BL+ S(6+ Lk), w= 2HL+ 6+ Lk. Then 
(A) becomes a quintic in k, the coefficient of k® being — L’. 

Third, if n = L = 0, (F) and (E£) give mandec. If m + 0, (C), (B), (A) 
give k, j, h, so that there is an unique rational solution. 

8. There remains the case b = g = 0, a +0. Multiplying x by 1/a, 
and y by a?, we may set a= 1. Hence we examine conditions (14) when 
b=g=0,a=1. 

First, let J + 0. Replacing w by w — zH/(27), we have H = 0. Then 
equations (14), other than the first two, give 


L = 0, p= n=—e, d= —G, k= (?+ F)/J, 
c=je—-@—-D, h= (C+ 2ce+ km)/J. 
We retain the abbreviations / and p = 2Ge — E+ ke, whence m = Jj + p. 
Equations (B) and (A) of (14) become 

JP —Ej+ta=0, a= B+ 26(4+4+ D) — 2ce*)/J — kep/J, 
F)P+ Bit+y=0, B= 2%2k+ C+ — D), 
y= At (C+ (CH 2ce)p/J + kp*/J. 
If in each coefficient of these two quadratics in 7 we retain only the term of 
maximum degree in e, we find that their resultant becomes 5e'*/J®. Hence 


e is a root of.an equation of degree 16. 
Second, let J = 0. Then equations (14), other than the first three, give 


L = 0, p = 0, n= —e-—H, d= — G, ec’ = — F, 
c=j(e+ A) m= k(e+ H) + 2Ge — E. 


The possibility of a (unique) solution with e = — H is easily decided since 
equations (A), (B), (C) determine uniquely h, 7, k if m+0. When 
e + — H, e + 0, we express the unknowns in terms of & and e, retaining 
the abbreviation m. By (C), c = (— C — km)/(2e). Equating this to 
the above value of c, we get 7. Then (B) gives h. Hence (A) becomes 


4Ae*(e + H)? + (C+ km)(e + H)? + 4Bme?(e + A) 
+ 4(C + km)Gme(e + H) + 2me{2(G? + D)e — C — km} = 0, 


which is a quartic in k, the coefficient of k* being (H + e)*(H — e). 

9. Consider the case, excluded in § 5, in which /3, /g, lz, ls are all free of z 
and hence are multiples of w. Postponing to the end of § 11 the case in 
which (11) is the product of y by its minor, we may assume, in view of the 


{ 
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substitutions S and T of § 3, that /; is the product of w by a constant, not 
zero, Which may be removed as a factor from the last row and multiplied 


into the last column of (2). Thus /; = w. We may take /; = 0 by (4).. 


Then — /3w = az?+ bw? by (3). Thus a= 0 and hence 6 = 0 by $5. 
Then /;= 0. Then /, = 0 and by removing its coefficient from the second 
row and multiplying it into the second column, we may set ; = w. Hence 


f= ls x = lL ow | = ay — + bls) + 
w 0 y | 


For y = 0, f reduces to wd, if 2 = hy +. Hence in 
fs=—yQ+ Kk, Q = Ay? + 2Byz + 2Cyw + Dz? + 2Ezw + Fu’, 


the cubic function K of z, w must have the factor w”, the quotient being 
X= jz+ kw. Further, Q must be of the form [| + lol, — kw?. We shall 
examine the last question independently of our main problem. If j + 0, 
we introduce the known function \ as a new variable z; let @ become 
Q’ = A’y?>+ ---. By (5) we may subtract a multiple of /, = hy + z from 
1, and assume that /; = cy + ew lacks z. When Q’ is divided by /2 to give 
a remainder free of z, the quotient gives /,; and the remainder must be 
[i — hw. Hence the latter must equal the value of Q’ for z = — hy, the 
conditions for which are 


c? = A’ — 2B’/h+ ce = C’ — eP=ht+F’. 
Hence the cubic 
(A’ — 2B/h + D'h)(h+ F’) = (C’ — E’h)? 
must have a rational root h such that h + F’ is a rational square e?. 
Next, let 7 = 0. Then k + 0 since we assume that y is not a factor of f. 
We may take'l, = cy + dz free of w by (15). When Q and /j — hw” are 
divided by /,, the remainders must be equal and a comparison of the quo- 


tients determines /,. Thus the function obtained from Q by replacing w 
by — yh/k must be identical with (cy + dz)? — ?h'/k?, whence 


h? h 
P=D, d=B-—“, &=A-2° 


10. Consider the case v = 0, excluded in § 5, in which /g = 0, 77 = z+ ww. 
By (3), 3 = — a(z — uw), b = — aw’. 
If a = b = 0, determinant (2) now equals 


(36) + Isls) + 
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This shall equal f = — yQ+ K, for Qand K asin§9. Thus K = 
We postpone to § 11 the simple case K = 0. One of the factors of K deter- 
mines /; = z+ uw. The other two factors determine \ and /, up to constant 
factors; but our determinant (36) is not altered when /, and J, are multiplied 
by t + 0 and I, = hy + X is multiplied by 1/t. Hence we may regard X, 
J; and J; as fully determined. We introduce \ as a new variable z and 
proceed as in § 9. 

Next, let a + 0. In view of (5) we may take = rw. If r + 0, we 
apply the substitution (2/4) (/3l7) which corresponds to the interchange 
of rows and columns, and have the case v + 0 treated in § 5. Ifr = 0, the 
determinant is identical with (1) if . 


— + bls) — ali(z + uw)(z — uw) = fs, 


a condition which is exactly of the type last discussed. 

11. Consider factorable cubic forms yQ, hitherto excluded. If Q itself 
has the factor y, yQ equals a determinant whose elements outside the 
diagonal are all zero. In the contrary case, we may take Q = ca?+ ---, 
where c + 0, after applying a linear transformation on 2, z, w with rational 
coefficients. Replacing cy by a new y, we have yQ, the coefficient of 2? in 
( being unity. Making a suitable addition to z, we obtain Q = 27+ q, 
where q is a quadratic form in y, 2, w. Thus yQ is of the form (1) with 
fe = 0, fs = 

First, consider equations (14) when a=bD=G=H=J=L=0, 
g +0. The last seven give at once 


p=j=0, k= 2d, n= —2e, gm=F+e’, gh= —D—@, 2gc= E — 2de. 
Multiplying the first equation (14) by — g’ and the next two by g, we get 
dE + 2eD+ gB=0, 2dF+eE+ 9C=0, D+ dek = Ag’ + fA, 


where A = LE? — 4DF. Multiply them by e, d, — 2, respectively, and add. 
We get 
2g7A+ geB+ gdC = — 3A, 


gB+ 2g9eD+ gdk = 0, 
gC + + 2gdF = 0, 


the last two being our first two equations multiplied by g. The determinant 
of the coefficients of g, ge, gd is 8 times the determinant 6 of 


(37) q = Ay? + Byz + Cyw+ D2? + Ezw + Fu’. 
Hence 86g? = — 3A(— A). Thus if 6 +0, A +0, the equations have 


solutions g, e, d with g + 0, and g will be rational if 6isa square. If A = 0, 


| 

t 
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the equations have no solutions with g + 0 if any two-rowed minor is not 
zero. But if those minors are all zero, there is a solution with g + 0 
except when g = Ay’, A + 0. Hence the cubic y(a?+ q) has a rational 
determinantal representation of the type in § 5 only when the determinant 
of g is a rational square + 0 and its minor EL? — 4DF is not zero, or in the 
trivial case when q is a perfect square. If we make the substitution 
z=Z+eY,w= W+dY, y= gY, we see that gq becomes DZ? + EZW 
+ FW? — 4AY? and that our determinant of type (11) reduces to the 
product of g by (38), which is the value of the initial determinant when 
d=e=B=C=0,g=1. 

Milder restrictions are imposed by the method of § 10 for a = b = 0, 
K=0. Thus vA/;= 0. If ls = 0, the determinant equals the product of 
y by its minor, a case treated below. Hence A = 0, l2 = hy. The case 
h = 0 is of the type just postponed. Writing hl, as a new J; in (36), we 
may take h = 1. In view of (5) we may assume that /; is free of y. Our 
cubic is of the form (36) if ¢q=—Jli— yl4+ Al; Hence, by (37), 
—1,= Ay+ Bz+ Cw. It remains to choose 


lL, = dz+ ew, ls = rz + sw, l, = vz + uw 

so that Dz? ++ Ezw + Fw? = — li + Al; The conditions are 
wr = 2+ D, ur + vs = 2de+ us= e+ F, 
These are linear equations in r, s, — 1, the determinant of whose coefficients is 
u?(d? ++ D) — uv(2de + FE) + v(e? + F) = 0, 
—wD+ wk — = (ud — ve)". 

The problem is solvable if we can choose rational numbers wu, v, not both 
zero, such that the left member* is the square p? of a rational number. 
Then rational values of d, e may be chosen so that ud — ve = p. Hence 
y(x? + q) has a rational determinantal representation of the present type 


if and only if 2 + q vanishes at a rational point having y = 0. 
If our cubic is of the form in § 9, we must have \ = 0, and 


q= — — + hw’, 


where we may assume that h + 0 and that /; is free of y (as in the preceding 
case). Thus — hyls equals the sum of the first three terms of (37) and 
— [i+ hw’ equals the sum of the last three. The conditions on /; = dz-+ ew 
are d? = — D, 2de= — E,h =e? + F. These determine d, e, h rationally 
if — D is a rational square + 0, and are satisfied if D = EL = 0. In the 


* That this condition is necessary is seen by taking z = u, w = — v in the proposed 
identity, whence /; = 0. 


| 
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respective cases, x? + q vanishes at a rational point having y = w = 0 or 
x=y=w= 0. 

It remains to consider the frequently postponed case of the rational 
representation of a product of a linear form / and a quadratic form Q as a 
determinant in which the elements of a row or column are 0, 0,/. Thus Q 
is to be represented as a two-rowed determinant. An evident necessary 
condition is that Q vanish at a rational point. Then, as in § 2, Q can be 
transformed rationally into zy + q(z, w) or q(y, 2, w). Any two-rowed 
determinant equal to the former may be given the form 


Cy!’ 
where B and C are free of x and y, whence A = 0, g = — BC. Similarly, 


if g(y, 2, w) equals a rational determinant, it vanishes for rational values 
not all zero of y, 2, w, and is rationally equivalent to yz + kw” or to a binary 
form. In the following summary, the various cases are presented in reverse 
order. 

THEOREM. A product of a linear form | and quadratic form Q in four 
variables with rational coefficients can be expressed as a determinant whose 
elements are linear functions with rational coefficients only in the following 
cases: (i) Q is expressible rationally as a two-rowed determinant vf it represents 
a cone which vanishes at a rational point not the vertex, or of Q vanishes at a 
rational point P and the tangent plane at P cuts the surface in rational lines. 
(ii) Q and 1 both vanish at the same rational point. (ili) 1Q is rationally 
equivalent to y(x? + q), where q is the ternary form (37) whose determinant is 
the square of a rational number + 0 and its minor A = E? — 4DF is not zero; 
the resulting representation is derived from the special case 

a+sEY —DY W | 
(38) Y(2?+ DZ?+ EZW+ FW? — 4AY?) = FY -—Z | 
Z W 

by a linear substitution of typez = Z+eY,w= W+aY,y= gl. 

Coroutuary. If A, D, F are rational numbers for which 

Q=2?+ AY?+ DZ? + FW? 
vanishes at no rational point, YQ is representable rationally as a deter- 
minant if and only if ADF is a rational square. 

I have found a simple proof of this Corollary independently of the 
present theory. 

12. In view of §§ 5, 9, 10, 2’y + Gz? + LW? is expressible rationally in 
determinantal form only when Gz* + LW? is a product of three factors with 
rational coefficients whence GL = 0. But if K is any homogeneous cubic, 
x’y — K(2a, 2, w) = 0 evidently has the solutions 


| 
| | 
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pA’, pK(A, B, C), pA*B, pAc, 


and no further rational solutions when x + 0. For, if A is any rational 
number + 0, we may define p by x = pA’, B by z = pA°B, C by w = pA°C; 
then y has the value specified. 

These solutions do not imply a determinantal representation of the 
surface since there do not exist three linearly independent linear relations 
between x, ---, w with coefficients linear in A, B, C, provided K does not 
have the factor x. In fact, all ‘such relations are linear combinations of 
Be — Az = 0,Cx — Aw = 0. But if K = 7’2, for example, we may secure 
such relations by taking new parameters A, B, D = A?/C; then, for 
oD/A’, 

z= cAD, y = oB’, z= oBD, w= 
Bz — Az = 0, — Ar+ Dw = 0, — Bz+ Dy = 0. 

13. The question whether the special cases postponed in § 5 are really 
exceptional is best answered by a different approach to the problem of 
normalizing our initial determinant. Its matrix is M=xA+yB+2C+wD, 
where, by (2), 


10 0 h & 0). 
A — 0) 1 0 B ts ty 0 
0 0 0 0 0 1 


t; being the coefficient of y in/;._ To further normalize M, we have available 
any constant matrices P and Q, whose determinants are not zero, such that 
in PMQ the coefficient of x is our A, while that of y is a matrix B’ of type B 
with t; replaced by t;. Thus PAQ = A, PBQ = B’. Set R= PP. Then 
AQ = RA, BQ = RB’, which are easily seen to require that 


[a B 0] 
R=Q=\7 6 Q|. 
0 p 


Thus PMQ = RMR, so that the normalization of M must arise by trans- 


formation by a matrix R of the form just given. Let ¢,; and d;; denote the 
elements of C and D in the 7th row and jth column. By (2), ¢33 = d33 = 0, 
C22 = — Ci, dee = — du. 

The first case postponed in § 5 is that in which ¢3, ¢23, ¢31, and C32 are 
all zero. Since R7'!CR then has the same four elements zero, this case is 
truly exceptional. Excluding it, we may take* c3; + 0. By choice of p, 

For 56 = 0, B = y = p = 1, transformation of C by R interchanges and 
C3, and c3»._ We also allow passing to the transposed matrix (with rows and columns inter- 


changed), thus treating one of two similar problems. Note that the transposed matrix 
cannot be obtained from C by transforming by a matrix of the specal form R. 


| 
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we may take c3;= 1. Transforming by R with a=6=p=1, y=0, 
= — C32, we may now have c3. = 0 as well as c3i1 = 1, ¢33 = 0. If also 
the last row of R“CR is 1, 0, 0, we must have 8 = 0,a= pp. The last row 
of R“DR is then 


6 
+ - dso, 0. 
p p 


Hence the special case in which dz. (denoted by »v in § 5) is zero is truly 
exceptional. Excluding it, we may choose y/p and 6/p so that d3; = 0, 
djg = 1. Then if R is such that the last row of R“DR is still 0, 1, 0, we 
have y = 0, 6 =, so that R is a similarity-matrix having the diagonal 
elements equal and having the remaining elements all zero. Thus R 
transforms every matrix into itself and no further normalization of B, C, D 
is possible. 

14. The difficult part of our problem is to identify the final determinant 
in (11) with any given ternary form (13) in y, z, w. The matrix of that 
determinant is M = yY + 2Z + wW, where 


d j [e k g 
Y=|m —e Z=|n —e —b 
0 oO 1 1 O 0 lo 1 0 


Let Y’, Z’, W’ denote the similar matrices in c’, ---, g’, but with the same a 
and same b. Do there exist matrices P and Q with constant elements of 
determinants not zero such that PMQ = M’? If so, the determinants of 
M and M’ differ only by a constant factor which we may assume is unity. 
It is more convenient to treat WQ= RM’, where R= P-. There are 
really several questions, depending upon what is assumed to be given. 

First, let only the determinant of M be given and require all matrices 
M, M’, --- of our special form and investigate their equivalence. Since 
this is our initial difficult problem with a supplement, we will illustrate the 
facts by means of an instructive example. Let |M| = w*— 23 — yew, 
so that E = G = — 1, L = 1, while the remaining coefficients of (13) are 
zero.* Thus (15)-(17) become 


— 82+ 6(3—g?) = 0, + €3— 9) =0, 5e(21— 97) + — 9’)? = 0. 


By the first two, = 0 or — (3—°)?/9. If = 0, then g = +1, 


=e=0. If +0, (2? — —9/4)=0. If g= +3, 6=—2, 
e=2. Ifg=+3/2, 6 = 1/4, € = — 1/4. The equations preceding (15) 
become 
gp = 9 = 1, gn = — 2e, gk = 26, gem = & + 26, 


gh = 2e — &, = 1 — — 2ée. 


* The only singular point of x?y + t; = 0 is (0, 1, 0, 0). 


i 
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Changing the sign of g merely changes the signs of the elements of the 
first two rows of M. Hence we take g = 1, 3, 3/2 in turn and obtain 


000 01 O 001) 
Y=Y’=|000!|, 
001 10 0] 010 

@ 

Of, 

1 0 0 1 


~ 


s(—yt2—w), ¢—yt4e+ Bw 

M” =| y+ 22+ 4w), — 

w y 

The conditions YQ = RY’, ZQ = RZ’, WQ = RW’ are satisfied if and only 
if 


1 —2 0] —1 2 
R=r|-2 1 d=r 2 Gi. 
2 0 8 


so that M and M’ are equivalent. But in WM” the coefficient Y” of y is 
of rank 2, while Y is of rank i, so that M and M” are not equivalent. 
Hence two matrices of type M with the same determinant may or may 
not be equivalent. 

A second question relates to a possible simplification of our initial 
problem in advance of its solution. Let a and b be given, while the remain- 
ing parameters c, h, ---, g in M are indeterminates. Can we find matrices 
R and Q with elements independent of c, ---, g such that MQ = RM’ for 
suitably determined elements c’, ---, g’, a, b of M’? If we can find such 
matrices R and Q not both similarity matrices rJ and gI, we can employ 
them to normalize M formally in advance of its computation. Denote the 
elements of the 7th row and jth column of R and Q by r;; and q;; respectively. 
By the third columns of YQ = RY’, 


eqis + = 113, mqiz — = 123, 


identically in c, h, m. Hence qiz3 = q23 = 113 = 123 = 0. By the third 
elements of the third rows of ZQ = RZ’ and WQ = RW’, 


By hypothesis, their determinant ab + is not zero. Hence = = 0. 
Then by the third rows of our matrix products, 


Gir = = = 133, +9), Q = 


| 
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Since we may replace R by r33'R, we may take r33 = 1. We now employ 
matrix P = R7 = (p,;), which has p13 = po3 = 0, p33 = 1. Then PM=M"’. 
In 


mpi hpi — cpr 0 
PY = | cpa t+ mpoo hpe— cp 0| = Y’, 
Cp31 + hpsi — cpse 1 
we see by the last row and the sum of the diagonal elements that 
Then pi; = 1 by |M| = |M’|. For pu = +1, we have the trivial case 
P=Q=I. For py = — 1, the multiplication of M by P on the left 
is equivalent to changing the signs of the elements of the first two rows. 
But the elements a and b in Z and W are to remain unchanged. Hence 
formal normalization of M is possible only when a = b = 0 and then consists 
in changing the signs of the elements of the first two rows. Hence the equation 
for g then involves only even powers. 


| | 
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THE IMPOSSIBILITY OF EINSTEIN FIELDS IMMERSED IN 
FLAT SPACE OF FIVE DIMENSIONS. 


By Epwarp KASNER. 


By the theory of quadratic differential forms it follows that a general 
riemannien manifold of m dimensions can always be regarded as immersed 
in some flat space of n dimensions, where n does not exceed 3m(m + 1). 
Thus if m = 4, as in the Einstein theory, the form 

ds? = Dgindxidxy, (1, k = 1, 2, 3, 4) 
can be immersed in an n-flat, where the possible values of n are 4, 5, 6, 7, 
8, 9, 10. | 

If now we require the manifold to obey Einstein’s equations of gravita- 
tion, G;, = 0, the question arises which of these values of n are actually 
realizable. The case n = 4 is trivial, since then the curvature vanishes 
and there is no permanent gravitation. We wish to show now that the case 
n = 51s impossible; that is, no Einstein manifold can be regarded as imbedded 
in a five-flat. 

In a flat space of five dimensions let rectangular coordinates be denoted 
by 21, %2, 73, 24, w. Then any four-dimensional manifold M in that space 
may be defined by a single finite equation 
(1) w = X3, 

The element of length of M is 
ds? = dai + + dx} + dzxi + dw’, 
where 


dw = fidzi, (J = 


is the total differential, summation with respect to the repeated index 1 
being understood as usual. This can be written 


(2) ds? = 
where 
(3) gi =1t+ fi, giz = fife, etc. 


The determinant of (2) is 


i fife ifs 
P | fife fofs fof 
| fifs fofs fafs 1+fi 
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which reduces exactly to 
(4) g=1t+fitfitfs + fi. 
The minors of (4’) also simplify, giving 


pa 


g 
(5) 


Our problem is to find the function f so that the ten functions (3) shall 
satisfy the ten gravitational equations* G;, = 0. 


CALCULATION OF THE CHRISTOFFEL SYMBOLS. 


We now calculate the Christoffel symbols (of the second kind) 


(6) {aB, }= (Jae, B + IBe, a «)- 
From (3) we find the derivatives 

gik, = (Sif) = fives + 
this formula holding whether the subscripts are unequal or equal. Sub- 
stituting in (6), we find 
Here the second derivative f,, appears as a factor, and it remains to carry 
out the summation f,.g”* involving only first derivatives of f. If the vari- 
able index ¢ takes the value y, we have 


(no summation on left), 


=f, 


where the prime indicates that e’ is to vary through the three values different 
from y. For the three values of ¢ different from , we have 


e’ Je! 
Therefore 
22'f, fe 
fg 
Hence the final value of our symbol is 
faa 
(7) {aB, vy} =" 
where 
(7') Py 1+ 22 = 2g — 1 — 


* For the notation see Weyl or Eddington, or § 1 of our paper “Einstein’s theory of 
gravitation: determination of the field by light signals,” this JouRNAL, vol. 43 (1921), p. 20. 
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CALCULATION OF THE EINSTEIN TENSOR G;;,. 


Take first the case of unlike indices, 
(8) Giz = {la, B}{28, a} — {12, a} L, — {12, a}. + Lr, 


where L = } logg. We now simplify the rest of the discussion by taking 
our axes in the original five-flat so the w-axis is normal to the four-spread (1) 
and so that the axes 21, 22, %3, 24 agree with the principal directions of the 
spread at the selected point, which we select as origin. This means that 
the partial derivatives of first order fi, fo, fz, fs all vanish, and that 


fu=fis = fi = fos = fu = Ju = 9. 
The four-spread (1) may therefore be written in power series form 
(9) w= + + + + higher terms. 
Thus at the origin 


these values defining the principal curvatures. We have also g = 1, so 
that L = 0, at the origin. 

For this special choice of axes we see from (7) that all the Christoffel 
symbols vanish. This throws out the first and second terms of (8). The 
third and fourth terms involve second derivatives of f, but we find that both 
terms vanish. Hence Gj, and of course the other five combinations with 
unlike subscripts, vanish identically. 

Take now the case of like indices 


(10) Gir = {la, B}{1B, a} — ie {11, a}, + Li. 

Here the first and second terms vanish. The third term is found to be 
— ky(ki + ke + hs + hy). 

The last term reduces to , 

Hence* 


(11) Gy = — ky (ke + ks + 


Thus the conditions that the four-spread (9) shall be of the Einstein type 
(that is, obey Gi, = 0) are 


* From this we may find the known formula for scalar curvature 
Cf. Eddington, “Report,” p. 75, or Page, Trans. Connecticut Acad., Vol. 23 (1920), p. 408. 


We may also connect with the elegant interpretations of the tensor Gi, given by Herglotz, 
Leipziger Berichte, 1916 (not accessible), and Vermeil, Géttinger Berichte, 1917. 
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ky (ke + ks + ks) = 0, 

ko (ks + ky) = 0, 
(12) 

k3(ka + hit ke) = 0, 

kea(ki + ko + ks) = 0. 

On the other hand the conditions that (9) shall be flat (euclidean or 
homoloidal) are 
(13) hike = ky = kok = = = 0. 

The two sets (12) and (13) are easily seen to be equivalent, each meaning 
that at least three of the four quantities ky, k2, k3, k4 must vanish. 

This shows that in a flat space of five dimensions we cannot construct 
four-dimensional manifolds (other than the trivial euclidean type) which 
obey Einstein’s equations. That is, an Einstein spread representing a 
permanent gravitational field can never be regarded as immersed in a five-flat. 

In flat space of six dimensions, actual Einstein manifolds exist; in 


particular, the solar field discussed in the next paper. 


CoLumBIA UNIVERSITY, 
New YorK. 


| 


FINITE REPRESENTATION OF THE SOLAR GRAVITATIONAL 
FIELD IN FLAT SPACE OF SIX DIMENSIONS. 


By Epwarp KASNER. 


In the preceding paper it was shown that no solution of the Einstein 
equations G';, = 0 can represent a four-spread imbedded in a five-flat (except 
in the trivial case where there is no permanent gravitation, so that the four- 
spread is itself flat). It follows that the solar field can only be imbedded in 
a flat space of more than five dimensions. By the general theory of quad- 
ratic forms in four variables, the maximum number of dimensions required 
is ten. We shall now show that the requisite dimensionality for the solar 
field is actually six.* The finite equations, in six cartesian coordinates, of 
the curved four-spread representing the solar field are given. As seen in 
formulas (7) below, a hyperelliptic integral is involved. This spread may 
be described as a geometric model of the exact field in which we are living. 

The field may be taken in the usual Schwarzschild form 


(1) ds? = ydt? — y dr? — r°dé? — r* sin? 6d¢’, 
where 
2m 


here m is the mass of the sun, r, 6, ¢ are polar coordinates, and ¢ is the time. 
Introducing cartesian coordinates, x, y, z, we have the equivalent form 


2m r— 2m 
2 ds? = — dz’? — dy? — d22 — ——— dr? + — dé 
(2) r — 2m + r 
where 


In (1) we have four square terms all with variable coefficients. In (2) 
we have three squares with unit coefficients, and two squares with variable 
coefficients. Our object is to obtain an equivalent form with six squares, 
all with unit coefficients. Of course, if we admit imaginary transformations, 
all the coefficients can be made + 1; but if we keep to the real domain, a 
certain number will equal + 1 and a certain number — 1. 

The fourth coefficient in (2) is a function of r alone; therefore we intro- 
duce a new variable R so that 


* Other fields in six dimensions exist and will be discussed elsewhere. 
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2m 
dk = 
om! 
that is, 


(3) R= ~8m(r — 2m). 
Our form then becomes 


2 2 2 2 2 
(4) ds? = — dz dy” — dz? — dR? + Re ion dt”. 


The last term is of course not the square of an exact differential, and our 
general theorem on five-flats shows that we cannot obtain a form with five 
unit squares. 

Taking the last two terms of (4), which involve merely the two variables 
R and t, we can, by the usual theory of surfaces, replace their combination 


R? 
2 
(5) dS? = — dR? + 16m? —— dt 
by a sum of three unit squares. We readily find 
(6) dS? = dX? + dY? — dZ’, 
where 
y= _ Rsint t R cost 


y= 
VR2 + 16m?’ VR? + 16m?’ 


256m* 
fv (R? + 16m?) 


The form (4) may now be written 


(7) 


(8) ds? = — dz? — dy? — dz? + dS’; 


but of course it must be remembered that while dz, dy, dz are exact differ- 
entials, dS is not. In fact dS is the distance element of the auxiliary 
surface defined by (7), where X, Y, Z (or rather X, Y, 7Z) are cartesian 
coordinates in an auxiliary three-flat, and R and ¢ are gaussian coordinates 
on this surface. This transcendental surface (involving a hyperelliptic 
integral), which apparently has escaped notice, we shall designate by (S). 
It is obviously a surface of rotation. To really obtain a form with all 
exact differentials we substitute (6) in (8). This gives our 

THEOREM. The exact solar field may be written with six squares of exact 


differentials 
(9) ds? = — dz? — dy? — dz? + dX?+ dY? — dZ’, 


where X, Y, Z are defined by equations (7) in conjunction with (3) and (2’), so 
that they are known functions of the world coordinates x, y, 2, t. 
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Thus if we consider a flat space with six rectangular coordinates 2, y, z, 
AX, Y, Z (or rather, on account of the minus signs in (9), with coordinates* 
a, ty, 1, X, Y, iZ), the finite representation of the solar field is a curved 
manifold M, of four dimensions (hypersurface) situated in a six flat and 
defined by (7) with (3) and (2’). 

If on this M, we take the parameter ¢ equal to zero we obtain a sub- 
manifold M; which is peculiarly simple. From (7) we have X = 0, therefore 
the distance element of M; is, from (9), 


(10) — dz’ — dy? — d2?+ dY? — dZ2’, 
where 


R 256m* 
VR? + 16m? (R? + 16m?) 


This shows that M; is in the five-flat (2, y, z, Y, Z). But going back to 
the form (4), and substituting dt = 0, we have as an equivalent element 
of M3, 

(11) — dx’? — dy? — dz? — dR’; 


or, changing all the signs, we use as the final form 
(12) do? = dx? + dy? + dz? + dR’, 


where R is defined by (3). This is obviously an M3; in a four-flat. It is 
easily seen to be a rotation manifold. In particular if we put z= 0, 
we have 


(13) dx? + dy? + dR’, 
which defines a paraboloid obtained by rotating the parabola 
(13’) R= v2m(r — 2m), 


supposed drawn in the (r, R) plane, about its directrix (the axis of R). 
This simple result, defining the geometry in a plane through the sun, was 
first obtained by Ftamm.f{ It illustrates the spatial measurements for a given 
value of the time parameter ¢. 
The surface (S) defined by (6) and (7) is obviously entirely different, 
and illustrates the relation of the space and time measurements. It is also a 
 *In Weyl’s terminology this would be described by saying that the flat space here 
employed has four negative and two positive dimensions. The six space is not properly 
euclidean, but what Hilbert (Géttinger Nachrichten, 1915, 1917) would term pseudo-eucli- 
dean. 
+L. Flamm, “Beitrige zur Einsteinschen Gravitationstheorie,” Physik. Zeitschr., 
Vol. 25 (1916), p. 168. A simple exposition is given in Weyl, ‘‘Raum, Zeit, Materie,”’ 
3d edition (1919), p. 228. 
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surface of rotation, but the generating curve is transcendental. From (7), 


we find 
where 
2 ae 
(14’) 


~F?+ 16m? 


If we take p as abscissa and Z as ordinate we have the required curve. 
Rotation about the axis of ordinates (with Z replacement by 7Z) produces 
the surface (S). The time parameter ¢ is represented by the angle of 
rotation. 

Flamm’s paraboloid, defined by (13'), and the new surface, defined by (14), 
thus supplement each other, and together characterize the metric of the solar 


field. 


CoLuMBIA UNIVERSITY, 
New York. 
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ON THE MOTION OF TWO SPHEROIDS IN AN INFINITE LIQUID 
ALONG THEIR COMMON AXIS OF REVOLUTION. 


By BIBHUTIBHUSAN 


INTRODUCTION. 


Though the problem of the motion of two spheres in an infinite liquid 
along the line joining their centers has been completely solved by various 
investigators, the first writer to attempt the corresponding problem for 
two spheroids or ellipsoids is Prof. Karl Pearson.* His method does not, 
however, admit of further development and therefore, does not lead to the 
complete solution of the problem. In a previous paper,} I have shown how 
the problem can be completely solved in the case of two spheroids of revolu- 
tion of small ellipticities, the motion of the solids being along their common 
axis of revolution. The present paper deals with the more general case of 
the same problem in as much as no limitation has been imposed as regards 
the ellipticities. It will be seen that the success of the problem depends on 
certain transformation theorems for spheroidal harmonics which were not 
known before, though the corresponding theorems for spherical harmonics 
were given long ago by Bessel. 

All the results in this paper are believed to be new. 


§ 1. Lemma. 


1. Let 01, and re, 62, w2 be two systems of spherical polar coérdinates 
having their origins at 0; and OQ, respectively and their polar axes in the 
directions of 0,02 and 020;; also let 21, yi, 21 and 2, yo, 22 be two parallel 
systems of rectangular cartesian coérdinates with the same origins, the polar 
axes being taken for the axes of x; and 22; let s denote the distance between 
O; and Oz. 

Since 

po = Vy3 + 22 = pi 


we have 


* “On the Motion of Spherical and Ellipsoidal Bodies in Fluid Media,’’ Part II (Quar- 
terly Journal of Pure and Applied Mathematics, Vol. XX). 
+ Bulletin of the Calcutta Mathematical Society, Vol. VII. 
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P,(cos 2) 1 dd 


(x2 + ips cos 


dd 

{8 — (a1 + cos 8) 
when |s| > |a1+ 7p: cos |, we can expand the right-hand side of this 
equation in a series of spherical harmonics,* so that 


P,,(cos (— 1)” 1 | (2m+1)P COs | dd, | 
ky dt” 


n! 


where t = s/k; and k; is any quantity. 

Following the notation of Sir W. D. Niven in his memoir “ On Ellipsoidal 
Harmonics ” (Phil. Trans., Vol. 128 (1891), A), let G(2) denote the external 
harmonics of a system of confocal spheroids having the center at 02 and 
0,0. as the axis of revolution. Then @(2) can be expressed in a series of 
spherical harmonics in the following way:tT 


Qnty! 0 0 D? 
(2n + iyi (=. OY2’ 2(2n+ 3) 


D* 1 
T 9-4. Qn + 3)(Qn + 5) 


where stands for (a3 — c3) (07/023), and dz, C2 (az > ¢2) are the semi-axes 


of any one of the quadrics of the second family. 


Now 
@ @ (2n)! Yrs 22) 4 
and 
H 22) = 12Pn(cos 2). 
Therefore 
(2n)! P,(cos 62) 
(Is) dd 


2°n! (Xe + cos #) 


Hence if k3 denotes a} — cz we get 


* Vide Todhunter, “The Functions of Laplace, Lame and Bessel,” p. 88. 

+ Vide Niven’s memoir, p. 245; Hobson, “On the Evaluation of a Certain Surface 
Integral and its Application to the Expansion in Series of the Potential of Elipsoids,”’ 
Proc. Lond. Math. Soc., Vol. XXIV, p. 91. 

¢ Niven’s memoir, p. 236. 
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2 1 
(2n + 1) + ips cos 


(n + 1)(n + 2) As 
2(2n+ 3) (x2 + cos 


__@n)! 1 
nl Jo Qn ke dd. 


Again 


2 ay 3 1 a+ cos 
G(2) (2n 1) n! ky ) as 


d” 
| 2(2n + 3) dint? | Qm(t). 


Therefore we have 


(2n + 1)! On (2 + tpe cos do 


2°*nin! Jo he 


LIAS ont yp, 


n+l 
n! kat) 0 m=0 ky 


a; dnt? 


d” 
where stands for 


We shall now express this in terms of spheroidal coérdinates. Supposing 
the spheroids are prolate, we write 


ky cos 6; cosh 


pi = ky sin 0, sinh = — — 1)!”; 


1 
=f 2 + cos ) ao 


if 0, (2 + cos ) = 
2 


0 


then 


and 


Substituting in (I) we get 


* 


= (2m + 1)1(m, (A) 
where 
2"n! d” gi 


If however the spheroids are of the oblate or planetary form, ai — ci is 

negative, so that we write y; = c; — a), and the corresponding codrdinates 
are 
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= COS 0; sinh 
pi = ¥1 Sin 4; cosh = yi(1 — + 1)”. 
Therefore (I) becomes 
= (— 1)™(2m-+ n)Pm(u2)Pm(tr2), (B) 


where 


Qn ! 12 n+2 
n +1 d 


(2n + 1)!7! dt’ 2(2n + 3) dtv*? | 


= Y2/y1 and t; = s/i71. 

(A) and (B) are the lemmas that will be needed in the discussion of the 
problems that we shall presently take up. 

2. We can put the expressions for ®,(m, n) and ®,(m, n) into more 
compact forms as definite integrals. For 


d*Qn(t) _(— Pu(p)dp 


Therefore 
+1 ar (n + 1)(n+ 2) ot 
t—p 
Hence 


Similarly we obtain 


@) (m, n) Qn Pate. 


§ 2. HyprRopDYNAMICAL APPLICATION. 


1. Let a1, (a; > and ag, (a2 > ¢2) be the semi-axes of two ovary 
spheroids having their centers at 0; and OQ, respectively, 0,02 being the 
direction of motion. 

Let be the velocity of 0; towards that of O2 towards also let 
s denote the distance between the centers 0,02 at any instant. Take du, 1, 
w and do, fe, w two systems of ovary spheroidal coérdinates having their 
origins at 01, O2 respectively, so that the surfaces of the given spheroids are 
= 10 and re 


* Todhunter, loc. cit. 


®,(m, n) = On Pn(p)dp. 
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The problem before us is to find a function ¢ satisfying the following 


conditions: 

(1) ; Vv’ = 0, throughout the liquid, 

(2) ¢@ = 0, at infinity 

(3) = — Pi(u1), when = Ajo, 
0 

(4) = — (u2), when = 


where ki = aj — ci and k} = a} — c}. 
A value of ¢ which satisfies the first two conditions is given by 


(5) = > [A nP + BrP n(u2)Qn(r2) 
n=1 
If by a proper choice of the arbitrary constants A,, B,, we can make it 
satisfy the other two conditions, then ¢ will be the required solution. 
By the Lemma (A), near the surface of 0;, we have 


= | + Be (2m + 1)@,(m, | 


n=1 m=0 
Substituting in the condition (3), we get 


= = | 


r= 


+ B (2m 0s) | 
where Q;(Aio) and P,(Aio) stand for the values of (0/0\;)Q,(\1) and 
(0/01) Pn(A1) when dio has been substituted for \;. This must be true at 
every point on the surface of the spheroid. Hence we can equate to zero 
the coefficients of the various zonal harmonics of 1. 
Equating the coefficients of Pi(u:) 
— = + 3 n) By. (6) 


n=1 


Equating the coefficients of P,(u1), p > 1, 
O = 10) + (2p + 1) &1(p, 0) Bn. (7) 


We shall have two similar equations from the remaining surface condi- 
tion which can, however, be written down from symmetry. Thus 
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— uke = B 191 ( (Aco) + 3 1, An, (8) 


n=] 


forp> 1 
O = 20) + (2p + 2(p, n) An, (9) 


n=] 


where ®,(m, n) can be written, by symmetry, from ®(m, n), viz. 


Thus we have four equations (6) --- (9) for the determination of the 


unknown constants A and B. 
2. The General Values of the Constants A and B.—Substituting the 


values of B; and B, (n > 1) from the equations (8) and (9) in the equation 


(6), we get 
1) 
— = AiQi Qi) — take + n)A, 


= 
— 32 ( (2m + 1)w1(1, m) 4 ,@2(m, n), 
or 
An (2m + 1)@,(1, m)@(m, n) 


@,(1, 1 


This equation can be written in the form 


OnmAn — 00) =c, (m= 1, 2,3, ad inf) (1) 


n=1 
where 


= (2m + m)@2(m, n)Pm(A20)/Qn (Azo), (2) 


and 
c= — 1)/Q:i (20). (3) 


In a similar manner, for the determination of B’s, we get the equations 


— 4BiQiQe0) = (m = 1, 2, 3, --- ad inf.) (4) 


where 
= (2m + 1)@(1, m)@(m, Aro); (5) 


= 1)/Qi (10). (6) 


Thus we have finally an infinite number of equations (1) and (4) for 
the determination of an infinite number of unknowns, A as well as B. 


and 
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The theory of the solution of the equations of this class has been worked 
out by Hill, Poincare, Koch, Toeplitz, Schmidt and others. Hence the 
values of the constants can be found so that the problems becomes deter- 
minate, the value of the potential function at any point being given by (5) 
Art. 1, § 2. 

The general values of A as well as B, so determined will, however, be of 
little help for numerical calculations which will be necessary so as to get a 
physical idea of the state of liquid at any point. So we give a method of 
approximating to their values. 

3. Approximate Values of the Constants A and B.—We have 


= (— n - — ——— 


Substituting the value of Q,,(¢) in series, viz. 


2(2m + 3) 


and rearranging the terms we get 


(m+n)! 
(2m + 1)! (2n + 1)! 


i (2m + 1)! (2n + 1)!2(2n + 8) gmtnts 
Thus the lowest order of ®;(m, n) is 


(= 


n) = 


(ki + 


central distance 


If the spheroids are so separated that we can neglect the terms of the order 
linear dimension 
central distance 


3 
) , we get from the equation (1), Art. 1, § 2 


Uk, 


We see generally that, for p > 1 
(p + 1)! 
Q,Qu0) 


Bp, = — A, + terms of higher orders. 
. . Pp 


A, 


by symmetry. 


A, = 


B, + terms of higher orders 


- | 
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We proceed to approximate to the values of the A’s and B’s 
First Approxvmation.— 


s* Bs, _ Q: 8 


Second Approximation.— 


2 1 


U1; 


and so on. 
4. Verification of the Results Previously Obtained. ae terms of 


the order (= eee | and higher, we have found that 
central distance 
vs > 
Expressed in terms of spherical harmonics 
Pi(cos 0; 
1 1 
3 1.2 Pa(cos 02) 
also 
212 + 
Qi(z) = 


Now if e; and e2 be the eccentricities of the generating ellipses of the two 
spheroids respectively, we know 
1/A1, = 1/Az0, ky = ke = 


Substituting in the equation (1) and neglecting powers of e; and e: higher 
than the second, we get as far as the required order 


| 
2 
ad 
| 
of 
a 
of 
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Now from the results of the previous paper, we have neglecting powers 
of e; and é2 higher than the second 


“61 
1 


Making allowance for the difference of notations between the two 
papers we have 


a= a;(1 — 


Therefore 
A; = — Zuai(1 — B,= 


linear dimension 


3 
Since the terms of the order (| ———— are neglected, 
central distance 


5 2 
Ay A3 Ag 


By 


ds}. \ ri 
Hence the equation (v) becomes up to the order of terms retained 


P3(cos 0 


1 


P, 0 
+ — $e2) — — + 


so that the results of the two papers agree at least up to t ler. Simi- 
larly it can be shown that the two results agree in their mos' general form. 

It may also be pointed out that, considered term rm, the series 
for the velocity potential function now obtained is more convergent than 
the series obtained in the previous paper. 

5. Conclusion.—In the preceding articles, I have studied the case when 
the two spheroids have the same common axis of revoliition which coincides 
also with the direction of motion. In a similar way, can be investigated 
the case when the axes of revolution of the two spliercids are at right angles 
to the direction of motion and the transformation theorem for spheroital 
harmonics that will be required can be obtained 

University or Caucutta, 1919. 
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